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ABSTRACT 

We present stellar kinematics for a sample of 10 early-type galaxies observed using the Space Telescope Imag- 
ing Spectrograph (STIS) aboard the Hubble Space Telescope, and the Modular Spectrograph on the MDM Obser- 
vatory 2.4-m telescope. These observations are a part of an ongoing program to understand the co-evolution of 
supermassive black holes and their host galaxies. Our spectral ranges include either the calcium triplet absorption 
lines at 8498, 8542, and 8662 A, or the Mg b absorption at 5 175 A. The lines are used to derive line-of-sight veloc- 
ity distributions (LOSVDs) of the stars using a Maximum Penalized Likelihood method. We use Gauss-Hermite 
polynomials to parameterize the LOSVDs and find predominantly negative /i4 values (boxy distributions) in the 
central regions of our galaxies. One galaxy, NGC 4697, has significantly positive central hA (high tail weight). The 
majority of galaxies have a central velocity dispersion excess in the STIS kinematics over ground-based velocity 
dispersions. The galaxies with the strongest rotational support, as quantified with vmax^stis, have the smallest 
dispersion excess at STIS resolution. 

The best-fitting, general, axisymmetric dynamical models (described in a companion paper) require black holes 
in all cases, with masses ranging from 10 6 5 to 10 9,3 M Q . We replot these updated masses on the - <r relation, 
and show that the fit to only these 10 galaxies has a slope consistent with the fits to larger samples. The greatest 
outlier is NGC 2778, a dwarf elliptical with relatively poorly constrained black hole mass. The two best candidates 
for pseudobulges, NGC 3384 and 7457, do not deviate significantly from the established relation between Mbh and 
a. Neither do the three galaxies which show the most evidence of a recent merger, NGC 3608, 4473, and 4697. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: kinematics and dynamics 



1. INTRODUCTION 

The wealth of new data from the Space Telescope Imag- 
ing Spectrograph (STIS) aboard the Hubble Space Telescope 
(HST) is dramatically improving our understanding of the cen- 
tral regions of early-type galaxies. Its long-slit design allows 
for more efficient measurement of galaxy kinematics than the 
previous Faint Object Spectrograph. The ~0"1 spatial resolu- 
tion of STIS resolves the influence of a 10 8 M black hole (BH) 
in an L* galaxy (velocity dispersion a = 200 km s" 1 ) out to ~ 
25Mpc. 

Our team has observed a sample of 10 nearby elliptical 
or lenticular galaxies using STIS in order to address many 
questions concerning the demographics of supermassive black 
holes. Are BHs present in the centers of all early-type galax- 
ies? With which properties of the host galaxy does the BH mass 
correlate? What can this tell us about the evolution of BHs and 
galaxies, and the relationship of inactive BHs in nearby galaxies 



to more distant active galactic nuclei? Our data can also address 
several issues in the evolution of galaxies. For example, why is 
there a dichotomy between power-law and core surface bright- 
ness profiles (Faber et al. 1997)? Were cores created by the 
scouring action of merging black holes? Why do a few ellip- 
ticals not fit well into either category (Rest et al. 2001)? Why 
are the properties of BH in bulges and pseudobulges so simi- 
lar, when these two types of stellar systems probably formed 
in quite different ways (Kormendy et al. 2002). Our sample 
contains galaxies with both core and power-law profiles. It 
also contains two possible pseudobulges: the low-luminosity 
SO NGC 7457 and the SO NGC 3384. 

The dataset presented in this paper has already been instru- 
mental in establishing a new, fundamental correlation between 
the BH mass and stellar velocity dispersion (Mm - a, Gebhardt 
et al. 2000 (GOO); Ferrarese & Merritt 2000). By combining our 
data with other published results, GOO demonstrate that the new 
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FIG. 1,— HST images of our 10 early-type galaxies. These are V-band (F555W) WFPC2 images for all except NGC 4564 (F702W), and NGC 4697 (WFPC, de- 
convolved). Each is a 200x200 pixel subimage, corresponding to 9" x9"on the WFPC2 and WFPC images. The label "C" indicates a core-type surface brightness 
profile (Lauer et al. 1995), while power-law galaxies are labeled only by their NGC number; NGC 4473 is intermediate between core and power-law. 



correlation has less scatter than the one between BH mass and 
bulge luminosity (see Kormendy & Richstone 1995; Magorrian 
et al. 1998; Richstone et al. 1998; and van der Marel & van den 
Bosch 1998). Here we give the details of data reduction lead- 
ing to this result. In a companion paper (Gebhardt et al. 2003), 
we discuss the modeling technique used to determine the BH 



masses, while Tremaine et al. (2002) discusses the slope of the 
Mi,h - o relation. 

The present paper complements the modeling paper by 
parameterizing the stellar line-of-sight velocity distributions 
(hereafter LOSVDs) with a Gauss-Hermite (GH) expansion. 
The parameterization of LOSVDs with GH polynomials has 
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become common practice over the last decade as increasing 
importance has been placed on the precise shape of the stel- 
lar LOSVD (van der Marel & Franx 1993; Bender et al. 
1994). High-quality galaxy spectra have demonstrated that stel- 
lar LOSVDs deviate significantly from a Gaussian distribution 
(Bender 1990; Gerhard 1993; van der Marel & Franx 1993). 
For edge-on, rapidly rotating galaxies, the LOSVD is often 
asymmetric. In other cases it has longer or shorter tails com- 
pared to a Gaussian. In the GH expansion, the first term mea- 
sures the mean velocity and velocity dispersion. The second 
term measures asymmetric deviations from a Gaussian, while 
the third term measures symmetric deviations from a Gaus- 
sian. Measuring the entire LOSVD, rather than just its first two 
moments (mean velocity and velocity dispersion), constrains 
the phase-space distribution of the stars and thus helps to re- 
duce the ambiguities in the mass distribution derived from stel- 
lar kinematics. Although the fully general LOSVD is used in 
our BH modeling, the GH parameters are a convenient means 
of describing an LOSVD. They provide a check on the orbital 
anisotropics determined by the modeling, and a means of com- 
parison to previously published results. 

The GH parameters can also be combined with photomet- 
ric and kinematic parameters to provide insight into elliptical 
galaxies. For example, Bender, Saglia & Gerhard (1994, here- 
after BSG) define mean parameters (h^) and (Ha) and find cor- 
relations between these and a^/a, v/a, M B , and v/a* for 44 
ellipticals. The greater generality of the GH parameterization 
has been exploited in deriving mass profiles from ground-based 
spectroscopy for NGC 4342 (Cretton & van den Bosch 1999), 
NGC 3115 (Emsellem et al. 1999), NGC 1399 (Saglia et al. 
2000), NGC 2974 (Cinzano & van der Marel 1994), and NGC 
2434 (Rix et al. 1997). In the case of M87, ground-based GH 
parameters have even been used to constrain the mass of the 
central BH (van der Marel 1994a). Many papers present GH 
parameters derived from ground-based spectroscopy for large 
samples of galaxies (e.g., Kaprolin & Zeilinger 2000; Kron- 
awitter et al. 2000; Fisher 1997). Checks and comparisons to 
our ground-based GH parameters come largely from the sam- 
ples of BSG and Halliday et al. (2001), which have 5 and 3 
galaxies in common, respectively. STIS now enables us to 
examine the GH parameters at sub-arcsecond resolution. For 
NGC 1023 (Bower et al. 2001) and M32 (Joseph et al. 2001), 
STIS has clearly revealed the influence of a BH on the stars. 

This paper proceeds as follows. We describe our galaxy sam- 
ple in §2, and the observations in §3. Section 4 details the re- 
duction of STIS and complementary ground-based data. In §5, 
we first describe how LOSVDs and GH parameters are derived 
from our spectra, and then we present our kinematics for indi- 
vidual galaxies. In §6 we will discuss our results. 

2. SAMPLE DESCRIPTION 

A collage of all 10 target galaxies is shown in Figure 1. 
The galaxies were chosen from those with reliable surface- 
brightness distributions at //irresolution. We attempted to in- 
clude a large range of luminosity and to exclude problematic 
objects such as those containing obscuring dust. (The galaxy 
NGC 4697 shows a dust disk which extends out to 3. "4, but the 
inner 1 ."0 does not appear to be seriously obscured.) 

The stepped grayscales are used in Figure 1 to show the va- 
riety of isophote shapes. Some galaxies, like NGC 821, have 
"disky" isophotes that will produce positive parameters a^/a, 
where «4 is the fourth cosine coefficient in the Fourier expan- 



sion of the radial deviations, and a is the semimajor axis of the 
isophote (Lauer 1985; Bender 1988). Others, like NGC 4291, 
have "boxy" isophotes that will produce negative 04/0. Boxy 
isophotes are typically found in galaxies with "core" surface 
brightness profiles (Faber et al. 1997). A "core" profile has a 
break between inner and outer slopes and the inner logarithmic 
slope, 7 = -dlogl/dlogr, must be less than 0.3. Similarly, 
"disky" isophotes are found in "power-law" galaxies, which 
have no significant break and have steeper slopes. Our sam- 
ple has six power-law and four core galaxies. We count NGC 
4473 as a core because of its surface brightness profile, but it 
shares many properties with power-law galaxies (see §5). The 
core galaxies are grouped on the bottom of Figure 1 . 

The most luminous core galaxies tend to have a large veloc- 
ity dispersion and low surface brightness, making measurement 
of the LOSVDs difficult. Thus, we selected three core galaxies 
with relatively high surface brightness (NGC 3608, 4291, and 
4473). However, we selected one example of a relatively low 
surface brightness core galaxy, NGC 4649, which has fi v ~ 
15.9 mag arcsec" 2 at O'.'l. It is comparable to M87 in luminos- 
ity, and it serves as a test case for measuring absorption-line 
kinematics in low surface brightness giants. 

3. OBSERVATIONS 
3.1. STIS Observations 

Table 2 gives specifications for the spectrograph/grating 
combinations used. Our STIS observations used only the 
G750M grating and the STIS CCD detector. The STIS CCD 
is a 1024x 1024 pixel CCD with readout noise around 3.8 e~ at 
a gain of 1.0. We binned the CCD by a factor of 2 along the 
dispersion axis in all of our observations to improve the signal- 
to-noise ratio (S/N). This raised our reciprocal dispersion to 1 . 1 
A pix" 1 . STIS has a spatial scale of 0."0507 pix" 1 in all of the 
configurations used here. 

Table 3 gives the details of the HST observations. Typically, 
6 orbits were devoted to each galaxy (« 2700 seconds per or- 
bit). However, the galaxies NGC 821, 3384, and 4697 had 11 
orbits, and the galaxy NGC 4649 had 18 orbits because of its 
low surface brightness. Each orbit was divided (CR-split) into 
two, 1350-s exposures except the first orbit of a visit, which had 
shorter exposures. For the first 2 visits of the program (NGC 
4473 and 821), the galaxy center was only dithered slightly (at 
the sub-pixel level) on the chip between orbits. Partway through 
the program we learned that the CCD had a rapidly varying pop- 
ulation of hot and warm pixels. We then requested wider sep- 
arations between dithers (20 pixels) so that a dark frame could 
be constructed from the data. Five dither positions were used. 
The second visit to NGC 821 used wide dithers, so only NGC 
4473 did not receive any wide-dither observations. 

For one galaxy, NGC 4697, two 1050 s exposures were taken 
with the G750M in setup 3 (Table 2). Here, the spectral range is 
centered near \ cen = 6581 A. We see Ha+[N II] emission lines 
which allow a useful comparison of gas and stellar kinematics. 
The reduction and analysis of these data is described elsewhere 
(Pinkney et al. 2003). 

Table 3 includes a K3 III and a G8 III star observed by us. 
We use these as LOSVD-fitting templates, for assessing aper- 
ture illumination corrections, and for measuring the STIS point 
spread function (PSF). We also used HR7615 from Bower et al. 
(2001). The spectrum of the G8 III star (HR6770) is shown in 
Figure 3. The template stars were observed through the same 
slits that we used for our galaxies. We used two slit apertures 



4 



in our STIS setups: 52"x0.1" and 52"x0.2". The aperture 
point-source throughputs are 64% and 77.8%, respectively, at 
8500 A. The O'.'l slit required a PEAKUP while the 0"2 slit did 
not. We used the 0."2 slit for the four core galaxies (NGC 3608, 
4291, 4473, and 4649) and the 0."1 slit for the six power-law 
galaxies. 

3.2. Ground-Based Observations 

We used the Modular Spectrograph (hereafter, Modspec) on 
the 2.4-m Hiltner Telescope at MDM Observatory for longslit 
spectroscopy (Table 4). The camera was equipped with either 
a thick, frontside illuminated LORAL 2048 2 CCD ("Wilbur"), 
or a thinned, backside illuminated SITE 1024 2 CCD ("Char- 
lotte"). The LORAL chip had 4.7 e~ readout noise, while the 
SITE had 5.45 e~. The Modspec Ca II setup produces compara- 
ble spectral resolution to our STIS + G750M setups (Table 2). 
Its spatial resolution is worse, of course, varying with seeing in 
the range 0"6-2"0. Our CCDs had pixel scales of 0."371 pix" 1 
(Wilbur) and 0"59 pix" 1 (Charlotte). 

Multiple exposures were taken of each galaxy, with lower 
surface brightness galaxies receiving more exposures. The slit 
was placed along the major axis and at least one other position 
angle (Table 4). Small dithers of the galaxy along the slit were 
used to remove CCD defects. The galaxy observations were 
bracketed by star observations to help monitor the seeing, and 
to build a library of template stars for defining the broadening 
functions in the galaxy spectra. Template exposures were gen- 
erally < 1 minute. Calibration frames included Ar comparison 
lamp spectra, continuum lamp spectra, and twilight sky spec- 
tra. Guiding used either starlight that was reflected off of the 
slit jaws or direct guiding from stars outside of the slit- viewing 
area. The root-mean-square (rms) deviations on the sky during 
the guiding were typically around 0.3". We often have a star in 
the slit along with the galaxy. We can compare the FWHM of 
the star during a long exposure with the FWHM of a star during 
a short exposure to check the guiding stability. For a 20-minute 
exposure, the guiding added about a 10-20% increase in the 
FWHM. 

4. DATA REDUCTION 

4.1. G750M Call Spectra 

The STIS data reduction was done with our own FORTRAN 
programs and FITSIO subroutines (Pence 1998). Our primary 
reasons for not using CALSTIS (Hodge et al. 1998) within 
IRAF are 1) a better dark frame can be created out of the data 
than the weekly darks used by CALSTIS, 2) many features of 
CALSTIS were not important to us, such as heliocentric cor- 
rection, 3) we did not want 2D rectification to be unnecessarily 
complex. 

We began by extracting the raw spectra from the multi- 
dimensional FITS file. To remove the bias level, we subtracted 
a constant fit to the overscan region. The overscan was about 
1110 data numbers (DN) for galaxy spectra, and 1505 DN for 
flats and comparison lamps, with little variation from one visit 
to the next. 

Subtraction of an accurate dark current was important be- 
cause the STIS CCD has warm and hot pixels which evolve on 
timescales of < 1 day. The mean dark current was about 0.007 
s -1 pix" 1 and this produced about 10 DN pix" 1 in our exposures. 
We used an archival dark only for the 2 visits that were not 
dithered (NGC 4473 and 821); for the rest, we created a dark 



out of the data. This self-dark required many iterations. The 
first dark was a biweight combination (see Beers et al. 1990) of 
the individual exposures after masking the pixels within ^25 
rows of the galaxy peak. The first dark, although imperfect, was 
subtracted from each spectrum. We proceeded with the flat- 
fielding, bad-pixel flagging, and un-dithering, until a 2D spec- 
trum can be made by combining the exposures. This 2D galaxy 
spectrum was then de-flattened, dithered, and subtracted from 
the raw spectra. The results contained the dark current, cosmic 
rays, and some residuals from the first, rough, dark subtraction. 
These were again biweight combined to form an improved dark. 
We then repeated reduction with this dark frame to form a sec- 
ond 2D spectrum which was again subtracted from the raws for 
a third self-dark iteration. At least 5 such iterations were per- 
formed to remove residual galaxy signal from the dark. A slice 
along the spatial axis showed the galaxy profile to dwindle to 
zero if low-level dark current was correctly subtracted from the 
final galaxy 2D spectrum. This dark also corrected for any bias 
pattern on the chip. The final galaxy spectrum was typically 
more free of deviant pixels than if an archival dark was used. 
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FIG. 2. — The ratio of two CCD flats compared with the original. The 
wavy line is an average of 6 rows across the middle of the combined flat for 
NGC 3384. It was taken with the 52"x0.1" aperture and the G750M 8561 A 
setup. The relatively flat line is the result of dividing this flat by another taken 
9 months later with the same setup. 

Within our self-dark iteration loop, the dark subtraction was 
followed by flat-fielding. The flat-fielding was important be- 
cause interference fringes become strong for A >7500 A. Our 
raw data had about 10% peak-to-peak variations due to fring- 
ing (Fig. 2). We used a flat which was a combination of the 
contemporaneous tungsten flats. These are intended to reduce 
fringe effects to the 0.9% level (Leitherer et al. 2001). It was 
difficult to measure fringe residuals in our galaxy spectra, but 
ratios of flats taken within a visit showed no detectable fringe 
pattern. Also, a ratio of flats taken 9 months apart through the 
52" x 0.1" aperture contained a residual fringe pattern at the 
<1.0% level (Fig. 2). The stability of these tungsten flats is 
one clue to their effectiveness in our reduction. Another is the 
shape of the final spectra (Figures 3 and 4), especially for stan- 
dard stars, which tend to show a linear continuum with only 
familiar features superimposed. The overall slope of the galaxy 
continuum does not have to be perfectly accurate for our pur- 
poses, so any color mismatch between the tungsten flat and our 
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FIG. 3. — Collage of spectra extracted from the central spatial bin for each power-law galaxy. For each galaxy, there is a pair of spectra: the top one is from 
Modspec and the bottom is from STIS. The spectra are the average (biweight) of all exposures. The STIS spectra have been normalized by a linear or parabolic fit to 
the continuum, while the Modspec spectra were divided by a fit to the local continuum in 100-pixel wide regions, with linear interpolation. The spectra are shifted 
vertically to avoid overlap. The template star HR6770 is shown on top for comparison. For NGC 821, there is no ground-based calcium triplet spectrum so we have 
shifted the Mg (5175 A) spectrum into the plotted wavelength range. 



galaxies is not problematic. 

After flat-fielding, the spectra had to be vertically shifted to a 
common dither, combined, and rotated. The shifting was done 
by measuring the peak of the galaxy profile in a 100-column 
wide region at the center of the CCD. The peak was measured 
with 1/9 pixel precision for each frame using a cubic spline 
interpolation. The spectra were shifted vertically to the same 
row of pixels, y=5 1 1 .0. Fractional pixel shifts were allowed 
in order to provide an accurate superposition, but this created 
smearing (redistribution of counts to neighboring pixels) along 
columns. The aligned spectra were then combined using the 
biweight to filter cosmic -ray hits. Since there were five dither 
positions, every final pixel value was a combination of at least 
five distinct pixels on the CCD chip. The resulting spectrum 
has a dispersion with a 0.6° clockwise tilt relative to the CCD 



rows. This was removed in the last step using a rotation about 
the chip center. The rotation adds additional smearing which is 
spatially periodic, i.e., regions of minimal smearing occur ev- 
ery ~95 pixels. The center of the spectrum, A=8559 A, was a 
smearing minimum which happens to fall near the Ca II 8542 
A line for our galaxies. 

The final 2D spectra contain imperfections. Figure 5 is an 
example of a final unsharp-masked, dithered dataset for NGC 
3384. The box overlay on NGC 3384 shows the ±1" extent of 
our ID extractions from the galaxy center. There is a low-level 
(c « 1 .0 DN) pattern noise apparent on Figure 5. It is possibly 
caused by residuals from dark-subtraction and/or flat-fielding 
which are then spatially repeated when the 5 dithers are com- 
bined. The dark contains charge trails along columns which 
become more numerous toward the bottom of the chip. These 
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FIG. 4. — Same as Figure 3, except here we show spectra from each core galaxy. 



features do not all subtract well leaving linear features on the 
spectrogram. Other imperfections are not easily visible. There 
are typically a dozen bright features between the fiducial bars 
for an average of ten, 1350-s exposures. These are 1-2 pixel 
events smeared by shifting to cover an area of <~ 4 pixels. They 
probably originate in pixels that contain cosmic rays in many 
of the exposures. (Cosmic ray hits occur on about 5% of pixels 
on the raw frames.) 

We extracted 1 -dimensional (ID) spectra from the final spec- 
trogram using a biweight combination of rows. Our standard- 
ized binning scheme (Table 6) used 1 -pixel wide bins near the 
galaxy center to optimize spatial resolution. All of the bins 
were contained within l."0 (20 pixels) of the galaxy center. 
Our dither pattern positioned the galaxy center on the rows 
y=471, 491, 511, 531 and 551, so that the residuals of the fidu- 
cial bars fell as near as 9. "2 to the galaxy center. These bars 
did not interfere with our spectroscopy because the S/N gener- 
ally becomes too poor to measure kinematics by r=l ."5 for our 



galaxies. Since our spectral extractions went out 20 pixels from 
the galaxy center and our maximum dithers were 40 pixels, all 
of our spectra originate between y=450 and 570 on the chip. 
This simplified several other aspects of the reduction, including 
wavelength calibration. 

A wavelength solution was obtained for the final galaxy spec- 
tra by summing the rows between y=450-550 on each of the 
Pt-Cr/Ne comparison lamp frames. There were typically six 
unique comparison spectra per visit. The final solution was an 
average of these six solutions. For each solution, we fitted a line 
to wavelength vs. pixel for 12 lamp lines and got residuals with 
rms«0.12 A. The lines were the 12 brightest after omitting a 
close double. A typical solution was y(A)=1.10896±0.00004x 
+ 8276.26±0.04, the errors are the standard deviation from the 
six solutions. Over the region occupied by the galaxy spec- 
tra (y=450-570), the lamp lines exhibit a slight bowing on the 
CCD which is asymmetric with respect to the central row. We 
considered errors caused by this bowing. The top dither at 



FIG . 5 . — Unsharp-masked 2D spectrum for NGC 3384. Only a subsection of the CCD is shown. The box is 2" tall and shows the vertical extent of the extracted 
ID spectra. The Ca II triplet absorption is visible in three vertical bands. 



y=551 has a zeropoint that differs by about +0.17 A from that 
of the central dither (y=511), and the bottom dither at y=471 
differs by about -0. 1 1 A. Thus, by combining vertically shifted 
dithers we artificially broaden the lines. The width of the lines 
increases by less than 0.5% as a result of this co-addition. Since 
the Ca absorption lines are broader than the lamp lines, their 
fractional change in width should be even less. The errors in the 
lamp-line centroids caused by the co-addition are small because 
the shift caused by the top dither will be nearly compensated by 
the shift from the bottom dither. Any remaining zeropoint error 
from the asymmetry is < 0.05 A. Finally, the galaxy spectrum 
was rotated by ~ 0.6° while the comparison spectrum was not. 
This introduced only negligible errors (<0.01 A) in line posi- 
tions. 

4.2. Template Spectra 

The template-star observations required a different reduction 
than the galaxies. Five exposures were taken through each of 
the 52"x0.2" and 52"x0.1" slits. The exposure times were 2.2 
s, so dark subtraction was unnecessary. Between exposures, the 



star was stepped across the slit by 0."04 and 0"025 for the 0"2 
and 0"1 slits, respectively. Combining each set of five spec- 
tra allowed us to simulate the observation of a diffuse source 
through both slits. The combined spectrum was then rotated by 
0.6°. Two lamp flats were obtained with exposure times of 142 
s (0."1 slit) and 84 s (0."2 slit). The flats were also combined and 
rotated. Next, ID spectra were extracted from both the star and 
flatfield spectrograms by averaging the 6 central rows. The ex- 
tracted flat was normalized and then divided into the extracted 
stellar spectrum. The resulting spectra (e.g., HR6770 in Fig. 3) 
appear to be corrected of fringing. 

The stepping of the stars across the slits allowed us to mea- 
sure the instrumental broadening of the Ca II triplet absorption 
lines. The first and last exposures in the 0"2 slit had the star po- 
sitioned at opposite extremes in the slit, a pointing difference of 
0."16 implying a 1.6 A centroid shift. However, the measured 
difference in centroids was only 1.07 A (38 km s" 1 ) in line cen- 
troids. The difference presumably originates in asymmetric vi- 
gnetting of the stellar PSF by the slit. The Ca lines had widths 
of about 3.4, 5.5 and 3.9 A, for 8498 A, 8542 A and 8662 A, 
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FIG. 6. — Galaxy light profile along the STIS slit (roughly major axis). Columns 200-300 on the STIS chip were averaged, corresponding roughly to the bandwidth 
8495 - 8605 A . Panel (a) shows the core galaxies. These were observed with the 52"x0.2" STIS slit. The "shoulder" in NGC 4473 is labelled (see §5.6). Also, 
dashed lines are overlayed showing the light in a simulated STIS slit over the WFPC2 data for NGC 4473. Panel (b) shows the power-law galaxies. These were 
observed with the 52" xO.l" STIS slit. 



respectively, in the individual 2.2-s exposures. The combined 
spectra showed only ~0. 1 A broadening when summed in the 
0"2 aperture. Thus, there is a < 5% broadening of the nar- 
rowest stellar Ca line due to the wider of the two slit widths. 
The broadening of Ca lines in our target galaxies will be even 
less. Finally, the STIS PSF was also derived from the template 
star exposures. The method and results are identical to those 
described by Bower et al. (2001). 

4.3. Ground-Based Spectra 

The Modspec data were also reduced with FITSIO routines. 
The CCD frames were overscan corrected, bias subtracted, and 
trimmed. The flat frames were constructed out of the twilight 
sky frames and the continuum-lamp frames as follows. We de- 
rived a small-scale structure flat from the continuum lamp. The 
large-scale structure along the dispersion axis was also taken 



from the continuum lamp. The structure along the spatial axis, 
however, was taken from the twilight frames. The final flat was 
the product of these three individual flats. 

After flat-fielding, the spectra were rectified. We used wave- 
length solutions from neon and argon comparison lamps and 
the traces of several standard-star spectra to define our geomet- 
ric transformations. Pixel smearing was not as much of a con- 
cern with the Modspec data as with STIS because we do not 
rely on Modspec for high spatial resolution. Therefore, nonlin- 
ear transformations were used and spectral images were shifted, 
aligned, and dispersion-corrected as needed before combining 
into final 2D spectra. Sky spectra were defined on the edges of 
the chip, where galaxy light was minimal, and subtracted from 
the 2D spectra. 

Table 4 shows our exposure times for each position angle 
(PA) for each galaxy. We took between 2 and 23 exposures for 
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a given PA to reach the desired S/N. As an example, a 1200-s 
exposure of an elliptical galaxy (NGC 4564 is used here) with 
an ffi>r-measured surface brightness at l."0 of /iy=16.0 mag 
arcsec" 2 will have a Ca spectrum with about S/N«20.0 per A 
for a spectrum extracted from a 0."37-wide bin (one pixel on 
the Wilbur CCD). For such a galaxy, 7 exposures will allow a 
S/N « 50 per A to be reached in the central extraction. This is 
easily sufficient for measuring the h3 and hA parameters of the 
LOSVD. The outer parts of the galaxy can be measured with 
adequate S/N by using broader extraction bins (see Table 6). In 
the case of NGC 4564, to reach S/N « 25 (giving uncertainties 
in H3 and H4 of ±0.05) at a radius where ^y=19 mag arcsec" 2 , 
a 4-pixel wide bin and 7 exposures are sufficient. 

4.4. Imaging 

A luminosity density distribution is needed for the modeling 
of each galaxy's mass distribution (see Gebhardt et al. 2003). 
A non-parametric fit is made to the surface brightness profile, 
and this is deprojected into luminosity density following Geb- 
hardt et al. (1996). We used primarily WFPC2 images in fil- 
ters F555W (V) and F814W (/) from HST proposals PID 5512, 
6099, 6587, and 6357. For NGC 4697, we used pre-COSTAR 
WFPC1 data in F555W. The images were typically sums of 4 
or more exposures. Starting with PID 6587, we used sub-pixel 
dithers for improved spatial sampling. All WFPC2 images were 
deconvolved with 40 iterations of Lucy-Richardson deconvolu- 
tion, while 80 iterations were applied to the WFPC images. Our 
images are shown out to r = 4. "5 in Figure 1 and surface bright- 
ness profiles are shown in Gebhardt et al. (2003). The profiles 
are actually a composite of HST and ground-based data, where 
the ground-based profile is shifted by a constant to match HST 
at a radius where seeing effects are minimal. The ground-based 
CCD photometry was taken primarily from Peletier (1989). 

For reference, the parameters of the best fit of the photome- 
try to the five-parameter Nuker law (Lauer et al. 1995) are given 
in Table 5. We show published values from Rest et al. (2001), 
Ravindranath et al. (2001), and Faber et al. (1997), or, if these 
are not available, our own fits to the / and V band WFPC2 im- 
ages. The Nuker law is a double power law given by 



/(r) = 2 (/3-7)/a 4 ^ 7 



1 + 



( 7 -/3)/a 



(1) 



where 7 is the slope as r — > 0, (3 is the slope at large r, and a 
determines the sharpness of the break between the two power 
laws at rb. These parameters are correlated (see Byun et al. 
1996), and the slippery nature of the 5-parameter fit allows for 
occasional sporadic values. These are fits to the major axis sur- 
face brightness distribution. If two sets of parameters are given 
for a galaxy, we adopt the average for each parameter for the 
purpose of plotting (§6). 

4.5. STIS Light Profiles 

Our STIS spectra also provide information on the near- 
infrared surface brightnesses of our galaxies. Figure 6 shows 
intensity as a function of slit position. The slits were approx- 
imately aligned with the major axes (Table 3). We produced 
intensity profiles in order to check the shifts of our galaxy 
peaks to a common position, row 5 1 1 .0. The profiles are slices 
along the spatial axis of the 2D spectrum, which are each a bi- 
weight of at least 10 exposures. The slices are 101 pixels wide 
(columns 200-300). The peak was defined as the maximum of 



a smooth, cubic spline fitted to the profile, thus allowing sub- 
pixel interpolation. The final 2D spectra were always peaked 
within 0.4 pixels of row 5 1 1 .0. For the creation of Figure 6, the 
fitted peaks were shifted to fall exactly on 5 1 1 .0 to best demon- 
strate any asymmetry in the galaxies. 

It is worth noting that these STIS datasets are valuable as 
probes of surface brightness features because total exposure 
times are over 10 times greater then typical WFPC2 surveys, 
the near-IR wavelengths are less dust sensitive, and saturation 
is avoided by dispersing the light. The ordinate of Figure 6 
gives the counts obtained in 1350 seconds ("DN" are equiva- 
lent to the number of electrons because gain = 1). However, at 
least 1010 pixels are averaged to determine each point so that 
the Poisson errors are actually very small: at 60 DN, the 68% 
error is only 0.2 DN, or about one-half the size of the plot- 
ted points. The kinematics are derived only out to ± 20 pixels 
(1 ."0). A point is labelled "double nucleus?" in reference to the 
note given in Table 1 of Byun et al. (1996). Here we find a sig- 
nificant asymmetry in the two sides of NGC 4473: one side is 
~ 5 DN higher than the other at \r\ =2 pixels. Also, the labelled 
point deviates from a line connecting the adjacent points by 1 
DN, which is also marginally significant. As can be seen in the 
steep profiles of the power-law galaxies (Figures 6b), the asym- 
metric PSF of STIS does not create asymmetries in these light 
profiles as strong as the one in NGC 4473. 

5. KINEMATICS FROM STARS: LOSVD FITTING 

The modeling of stellar kinematics requires the derivation of 
LOSVDs at numerous positions on the galaxy. The STIS spec- 
tra are used to sample only ±1"0 along the major axis, while 
our Modspec observations sample to 1-2 half-light radii (rang- 
ing from 10-70") along the major axis, and other position an- 
gles. After some experimentation, the binning shown in Table 6 
was chosen. It provides the best spatial resolution while main- 
taining adequate S/N. The first three bins are each 1 STIS pixel 
wide. 

Obtaining internal kinematic information requires a de- 
convolution of the observed galaxy spectrum using a tem- 
plate spectrum composed of several representative stellar spec- 
tra. Both the deconvolution process and the template li- 
brary are significant issues for obtaining the LOSVD. We de- 
convolve each spectrum using two different techniques: a 
maximum-penalized likelihood (MPL) estimate that obtains 
a non-parametric LOSVD, and a Fourier correlation quotient 
technique (FCQ; Bender 1990). MPL proceeds as follows. We 
choose an initial velocity profile in bins (the choice of initial 
profile has no effect on the result). The initial velocity profile 
convolved with a weighted-average template (discussed below) 
provides a galaxy spectrum that is used to calculate the resid- 
uals to the observed spectrum. The program fits the shape of 
the losvd by subdividing it into 13 velocity bins, and varying 
the counts in these bins while adjusting the template weights 
to provide the best match to each galaxy spectrum. The exact 
number of bins has little influence on the final results. Thus, the 
template mixture can vary from one radius to the next, just as 
the velocity, velocity dispersion, etc. can. The MPL technique 
is similar to that used by Saha & Williams (1994) and Merritt 
(1997). However, our method differs from theirs in that we fit 
simultaneously for the velocity profile and template weights. 

We continuum-divide all template and galaxy spectra. We 
use a local continuum for the estimate as opposed to a global 
n-parameter fit to the full spectrum. The continuum is esti- 
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FIG. 7. — A comparison of FCQ (circles) and MPL (triangles): kinematic profiles for NGC 4473. These parameters were derived from unsymmetrized spectra 
from STIS. From top to bottom, we show mean velocity, velocity dispersion, h3 (3rd coefficient from the Gauss-Hermite expansion), and h4 (4th coefficient). The 
solid lines are derived from the symmetrized LOSVD from MPL. 



mated by dividing the spectrum up into about 10 wavelength 
windows, estimating a robust mean (the biweight; see Beers, 
Flynn, & Gebhardt 1990) in each window, and interpolating 
between these values. This mean is measured from the high- 
est 1/3 of the points to compensate for absorption features. We 
have tried a variety of continuum estimates (varying the num- 
ber of local windows and number of points used in the local 
averaging) and find insignificant differences in the kinematic 
results. 

We use Monte Carlo simulations to measure the uncertainties 
on the velocity profile bins. For each realization, we generate 
a simulated galaxy spectrum based on the best-fit velocity pro- 
file and an estimate of the rms residual of the initial fit. The 
initial galaxy spectrum comes from the template star convolved 
with the measured LOSVD. This provides a galaxy spectrum 
with essentially zero noise (the noise in the template is insignif- 



icant for our purposes). From that initial galaxy spectrum, we 
then generate 100 realizations and determine the velocity pro- 
file, and hence the velocity dispersion, each time. Each real- 
ization contains flux values at each wavelength that are chosen 
from a Gaussian distribution, with the mean given by the ini- 
tial galaxy spectrum and the standard deviation given by the 
rms of the initial fit. The 100 realizations of the velocity pro- 
files provide a distribution of values from which we estimate the 
68% confidence bands. These velocity profiles and their 68% 
confidence bands are used directly in the dynamical modeling 
(Gebhardt et al. 2003). 

We wish to compare the MPL-derived LOSVD with those 
derived using FCQ (see above). To do this, we must convert 
the non-parametric LOSVDs into Gauss-Hermite polynomials. 
We use the least-squares estimator MRQMIN from Press et al. 
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FIG. 8. — Investigation of bias in MPL estimates of Gauss-Hermite parameters. From top to bottom, we show velocity dispersion, /13 (3rd coefficient from the 
Gauss-Hermite expansion), and A4 (4th coefficient). Each point is a different realization of a broadened template spectrum with noise added. The signal-to-noise of 
the input spectra increases to the right in each row of panels. In the top three panels, the colors light blue, blue, black, red and green represent /14s of -0.09, -0.04, 



0.0, 0.04 and 0.09, respectively. In the bottom six panels, green stands for a m 



-- 70 km s 1 , red for 170 km s 1 , and blue for 270 km s 



(1992) to find the best fit to this parameterization: 

f(y) = I exp(-y 2 /2)(l + h 3 H 3 (y) + h^iy)) , (2) 
where the mean velocity Vfu and velocity dispersion tr/„ are 
contained within y = (v- v/,>)/o/,>. Iq is the amplitude of the 
LOSVD at y = 0. The coefficient h 3 multiplies the asymmetric 
deviations from a Gaussian (i.e., skewness), while H4(y) pa- 
rameterizes symmetric deviations from a Gaussian (kurtosis). 
The definitions of Hj and H4 are given in BSG The first and 
second moments of the LOSVD will differ from v/a and a fa 



when hi and are non-zero. For example, the oy,, is ~ 10% 
smaller than the second moment about the mean when hi, = 0.1, 
and ~ 5% smaller for /13 = 0.1. Similarly, v/„ will be ~ 15% 
larger than the first moment for /13 = -0.1 (BSG). 

We tabulate the first four terms of the Gauss-Hermite expan- 
sion of the LOSVD for each galaxy in Tables 9 and 10. We use 
v/a and a fi t rather than corrected values (e.g., vo, c and <7o, c used 
in Joseph et al. 2001) for reasons listed in BSG, and for the most 
direct comparison to published results. Hereafter, we refer to 
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FIG. 9. — Kinematic profiles for NGC 821. From top to bottom, we show mean velocity, velocity dispersion, ft3 (3rd coefficient from the Gauss-Hermite expan- 
sion), and h4 (4th coefficient). Left: major axis data from r < I'.'O. The solid lines are derived from the symmetrized LOSVD from ST1S Ca II data. The circles are 
unsymmetrized data from STIS Ca II. The squares are ground-based, unsymmetrized data. Right: the entire radial extent of only the ground-based "Modspec" data. 
Again, a solid line connects the symmetrized, major-axis data, and the squares represent unsymmetrized, major axis (PA = 0°) data. Other symbols are labelled by 
the rotation of the slit with respect to the galaxy's major axis, in degrees. For plotting purposes, the prograde side is always on the right (r > 0). For NGC 821, both 
ground-based PAs are from spectra centered on Mg b. 



Vfu and a fit as v and a. We found that FCQ and MPL gave con- 
sistent profiles and errors. Figure 7 plots the first four Gauss- 
Hermite moments for both FCQ and MPL for NGC 4473. FCQ 
showed more scatter than MPL in galaxies with poorly defined 
lines (e.g., NGC 4649). We plot only the MPL-derived v, a, Kb, 
and h.4 against r in Figures 9-22. 

The Monte Carlo simulations used to measure the uncertain- 
ties also provides an estimate of the estimator biases. This bias 
may be important for the dynamical models since the shape of 
the velocity profile at some level determines the internal orbital 
structure. Furthermore, Joseph et al. (2001) argue that H\ es- 
timated from FCQ is biased to lower values when the galaxy 
dispersion is low (less than 100 km/s). The Monte Carlo sim- 
ulations provide the most accurate estimate of any biases since 



they simulate the exact instrumental setup and spectral sam- 
pling. We find little to no bias in our MPL estimates of the first 
four Gauss-Hermite moments. Another way to check biases 
and confidence bands is to fit to simulated datasets. We have 
done this for both FCQ and MPL. The FCQ results are similar 
to what is discussed in Joseph et al. (2001); that is, for galaxies 
with dispersion below 100 km/s, H4 is biased to lower values by 
about 0.04 for S/N=30 per angstrom, which is typical for these 
data. Figure 8 shows the results for MPL run on simulations 
using the same setup we used on the STIS galaxies, and for S/N 
ratios of 15, 30 and 60 per angstrom. The MPL results show 
no bias in either the Hi, or H4. Also, the uncertainties deter- 
mined from the simulated datasets agree with those determined 
from the galaxy spectra at the same S/N. I.e., at S/N=40 per 
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FIG. 10. — Kinematic profiles for NGC 2778 (see Fig. 9 for the meaning of the symbols). For NGC 2778, all ground-based PAs are from Ca spectra. 



angstrom, we obtain MPL uncertainties on Hi and H4 to around 
0.02. After correcting for the bias in FCQ, the uncertainties on 
Hi and H4 are similar. 

While we report the first four Gauss-Hermite moments in Ta- 
ble 10, we use the non-parametric velocity profile for the dy- 
namics. These moments are reported only to provide a sim- 
ple parameterization of the velocity profile for comparison and 
correlation studies. An important difference is that the non- 
parametric velocity profiles are not allowed to be negative at 
any point. Using only the first four moments of a Gauss- 
Hermite expansion will allow negative portions of the velocity 
profile. This is unphysical, but the effect should be small (at 
least when hi and /14 are small). Thus, the optimal compari- 
son for the dynamical models is to use non-negative velocity 
profiles. 

We now discuss each individual galaxy. We will include spe- 
cial points of interest, and a comparison to published kinemat- 
ics, where available. 



5.1. NGC 821 

NGC 821 is classified as . E . 6 . $ . (de Vaucouleurs et al. 
1991, hereafter, RC3). The uncertainty is likely to originate 
from its resemblance to an SO. NGC 821 is a clear example 
of an elliptical with disky isophotes (Lauer 1985; Bender et al. 
1988 (BDM); Nieto et al. 1991). As is typical of such galaxies, 
it has a power-law surface brightness profile at HST resolution 
(Ravindranath et al. 2001). 

NGC 821 is not detected in X-rays (Pellegrini 1999), or in 
Ha (Macchetto et al. 1996). It also does not appear to contain 
dust (Forbes 1991) although it is detected by IRAS in 100 ytm 
(Roberts et al. 1991). Trager et al. (2000) estimate the average 
age of its stellar population to be 8 Gyr. 

Kinematics show NGC 821 to be a fairly rapid rotator; Ni- 
eto et al. (1988, 1994) measure (v max la)* ~ 0.71 and we find 
~ 0.54 for our ground-based data (Table 8) and 0.75 using our 
higher STIS v max (Table 7). We get v max la ~ 0.4 (Tables 7 and 



14 



-1 -0.5 0.5 -20 -10 10 20 

=n i i i i i i i i i i i i i i i i i 1 1 i- -i J. i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i- 




-1 -0.5 0.5 -20 -10 10 20 

R (arcsec) R (arcsec) 



FIG. 11 .— Kinematic profiles for NGC 3384 (see Fig. 9 for the meaning of the symbols). For N3384, the dashed line (PA = 0°) and PA = 88° circles are Mg 
spectra. All others are Ca spectra. 



8) which is lower than previous values (0.49 in Bender 1988; 
0.65 in Busarello et al. 1992). 

Our kinematic profiles from Modspec are similar to those of 
BSG. Both show that the rotation rises quickly to ^60 km s" 1 
at 2". The STIS kinematics show the largest velocity gradient 
at 0."17 in the sample (Table 7). hi nicely mirrors the veloc- 
ity profile, as is typical in a rotating, edge-on elliptical. Our 
ground-based h4 shows a different behavior than in BSG: we 
find positive values at r > 2"and negative values in the cen- 
tral bin, whereas BSG finds only h4 < 0. Our values are prob- 
ably more reliable in this case because we have a longer ef- 
fective exposure time. The main difference between our data 
and others is the central rise in velocity dispersion. We mea- 
sure a higher ground-based value, 235±2 km s" 1 . McElroy's 
(1995) weighted average of published velocity dispersions was 
207 km s" 1 , and Hypercat (Prugniel et al. 1998) gives the aver- 
age of 12 velocity dispersions to be 209 km s" 1 . Our Modspec 



value does not exceed our central STIS measurement (248±16 
km s" 1 ), and may simply be a result of good seeing. When we 
measure the rms dispersion within a slit aperture of length 2r e 
(i.e., <7 e ), we find 209 km s" 1 , in agreement with the published 
averages. 

5.2. NGC 2778 

NGC 2778 is the second least luminous elliptical in the sam- 
ple with M B = -18.6 mag. It is located close to NGC 2779 in a 
group. Trager et al. (2000) estimate the ages of its stars to be 
only 6.0 Gyr. Its isophotes indicate an E2 morphology with a 
power-law profile (Peletier 1990; Lauer 2002, private commu- 
nication). 

Our STIS observations for this galaxy were the only ones 
taken in Backup Guiding Mode during the HST visit. During 
galaxy acquisition, only one of the two guidestars was found so 
that drifting about the roll axis was only constrained by gyros. 
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FIG. 12. — Kinematic profiles for NGC 3608 (see Fig. 9 for the meaning of the symbols). For NGC 3608, all ground-based PAs are from Ca spectra. 



Our parallel observations with the WFPC2 show virtually no 
change in the centroid of a galaxy image between the first and 
last (6th) orbit. Since the PC chip is about as far from the guide 
star as the STIS slit, there was probably no significant error in 
the STIS guiding. 

Davies et al. (1983), Gonzalez (1993) and Fisher et al. (1995) 
all provide kinematic profiles for NGC 2778. These rotation 
curves extend beyond ours to ^25 ."0. The maximum rotation 
velocity, ~ 100 km s" 1 , is reached by about 5". No rotation was 
observed along the minor axis; the galaxy is approximately an 
oblate, rapid rotator (Jedrzejewski & Schechter 1989). All dis- 
persion profiles show a steep rise from a « 100 to 200 km s"'in 
the inner 5". This galaxy has the greatest deviation from the 
Mhh - a relation in our sample: it has a small BH mass for its 
effective dispersion. The low surface brightness (see Fig. 6) and 
poorly resolved sphere of influence (GM^/er 2 ~ 0.02") make 
the BH detection in NGC 2778 relatively uncertain (Gebhardt 
et al. 2003). 



5.3. NGC 3384 

NGC 3384 is classified as SB (s) 0- : (RC3). It is in the 
Leo-I group, neighbored closely by the elliptical Ml 05 (NGC 
3379) and the spirals M95 and M96. It appears to have been 
named redundantly in the NGC as NGC 3371. Several H I 
clouds orbit around NGC 3384 and 3379 in a 200 kpc diam- 
eter ring formation (Schneider 1985). The H I has been inter- 
preted as gas stripped during an encounter between NGC 3384 
and the spiral NGC 3368 (Rood & Williams 1985). The inner 
regions of NGC 3384 appear largely free of interstellar matter 
(ISM): Tomita et al. (2000) find no dust in WFPC2 color-excess 
images. 

In the WFPC2 V and / images (Fig. 1), NGC 3384 appears 
consistent with a lenticular morphology. The Nuker-law profile 
shows a rapid transition (a = 1 1 .2) between inner cusp region 
(7 = 0.6) and the outer power-law (/3 = 1 .8) at the break radius 
rt, = 2."9 (Lauer 2002, private communication). There also ap- 



16 



-1 -0.5 0.5 -10 10 




-1 -0.5 0.5 -10 10 

R (arcsec) R (arcsec) 



FIG. 13. — Kinematic profiles for NGC 4291 (see Fig. 9 for the meaning of the symbols). For NGC 4291, all ground-based PAs are from Ca spectra. 



pears to be light in excess of the Nuker law at r < 0"3. Tomita 
et al. (2000) use the V and / images to search for dust and blue 
nuclei. They suspect a blue nucleus in NGC 3384, but claim 
that it is uncertain because the center is saturated on the PC. 
High resolution, ground-based imaging also hints at a nucleus 
(Kormendy 1985). Our deconvolved surface brightness profiles 
show a nearly fiat V-I profile. Our near-infrared light profile 
along the STIS slit (which is not saturated) shows a strong cen- 
tral peak, but not necessarily a stellar nucleus (Figure 6). Also, 
we do not see Paschen absorption in our Ca II spectra which 
would indicate a young stellar population. 

NGC 3384 has the highest S/N (> 40 in central bins) 
among our STIS datasets owing to its high surface brightness 
(A*v,o.i = 13.7 mag arcsec" 2 ) and 11 orbits of exposure (Table 3). 
Our ground-based kinematics are also of high quality (Fig. 1 1). 
We have both Ca II triplet (solid line) and Mg b (dashed line) 
data along the major axis. The results we obtain from these two 
wavelength ranges are consistent. 



Our Modspec kinematics for NGC 3384 can be compared to 
the integral field spectroscopy of de Zeeuw et al. (2002), and 
to long-slit spectroscopy by Fisher (1997), and Busarello et al. 

(1996) . All rotation velocity curves peak at r ~ 3."0 followed by 
a decrease. Then, beyond 10", the rotation climbs again. Such 
an "overshoot" is seen in other SO galaxies, most notably NGC 
4111, and in about 5-10 of the 18 SO galaxies studied by Fisher 

( 1 997) display some sort of overshoot or kink in their rotation at 
~ 5". De Zeeuw et al. and Fisher plot minor-axis profiles which 
indicate that the rapidly rotating disk component has a lower ve- 
locity dispersion than the bulge. This is also seen in our 75° and 
88° rotation curves which show a higher dispersion beyond 2" 
than the major axis data (Fig. 11). Thus, it appears that a cold 
stellar disk contributes differentially to the total light, increas- 
ing the fitted rotation velocity and decreasing the dispersion at 
small radii. At STIS resolution, however, the major axis disper- 
sion profile shows a distinct peak at r <0."2. The stellar disk 
also contributes to the strongly asymmetric LOSVDs: our data 
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show a peak h3 value of 0.15 along the major axis, the highest 
/i3 in our sample after NGC 4697. Our ground-based h4 val- 
ues are systematically lower than those of de Zeeuw et al. and 
Fisher, but all show a dip at r < 2". 

NGC 3384 has some of the observed characteristics of a 
pseudobulge, i. e., a high-density central mass component that 
resembles a bulge but that is thought to have been built by the 
inward transport of disk gas by (in the present case) a bar. 
This secular process is qualitatively different from the clas- 
sical picture in which bulges form on approximately a col- 
lapse timescale by violent relaxation and dissipation in a galaxy 
merger. Pseudobulge characteristics of the present galaxies are 
discussed in Kormendy et al. (2002), and more general reviews 
of pseudobulges are given in Kormendy (1993) and Kormendy 
& Gebhardt (2001). In essence, pseudobulges are left with 
some memory of their origin - they have disky structural prop- 
erties. In NGC 3384, the pseudobulge dominates the light at 



radius r < 6", i. e., interior to the bar, which is most impor- 
tant at r ~ 10" to 20" (Busarello et al. 1996; Jungwiert et al. 
1997). In NGC 3384, the evidence for a pseudobulge includes 
the highest v/er value in Table 7, an unusually high (v/er)* value 
in Table 8, and the observation that the "bulge" is almost as flat- 
tened as the disk (its ellipticity e ~ 0.45 in Table 8 as compared 
with e ~ 0.50 for the outer disk (Davoust et al. 1981; Busarello 
et al. 1996). Consistent with this, Bureau et al. (2002) and de 
Zeeuw et al. (2002) find evidence for a "disk embedded in the 
bulge" (high rotation and a flattened light distribution) from 
SAURON observations. 

5.4. NGC 3608 

NGC 3608 is a core galaxy in the Leo Group, separated by 
only 5.8' (39 kpc) from NGC 3607. It is undetected at 6 cm 
(Wrobel & Heeschen 1991) but contains extended X-ray emis- 
sion (log(L/ergs~ l ) = 40.4, Pelligrini 1999). Surface photome- 
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try shows an E2 elliptical with slightly boxy isophotes beyond 
5". Nuker-law parameter fits are published for the WFPC1 data 
in Lauer et al. (1995) and Faber et al. (1997), and we tabulate 
the parameters for the WFPC2 V image in Table 5. NGC 3608 
has a definite "core" profile with 7 «0.06. There is evidence 
for patchy dust in the WFPC2 data. Carollo et al. (1997) note 
an elongated dust feature in the innermost 0"6 (67 pc for d = 23 
Mpc) that looks like an off-centered ring. Tomita et al. (2000) 
estimate the dust feature to be 190 pc across. Singh et al. (1994) 
suggest that the interacting neighbor, NGC 3607, has acquired 
gas and dust from NGC 3608. 

Our ground-based kinematics can be compared to those in Je- 
drzejewski & Schechter (1988, 1989). Our mean velocity and 
velocity dispersion profiles are consistent with those published 
by these authors, although they plot a few more points beyond 
20". The outstanding feature of NGC 3608 is the kinematically 
distinct core. The slope of the rotation curve changes sign at 
about 5" (coincident with the onset of boxy isophotes), and the 
rotation itself changes sign at r m 13". The maximum rotation 
is only 24 km s" 1 , giving NGC 3608 the lowest v/cx in our sam- 
ple. Our STIS velocity dispersion profile shows a central peak 
which is higher than the ground-based value by about 100 km 
s" 1 . Also notable is the significantly negative central hA values 
in the ground-based data. 

5.5. NGC 4291 

NGC 4291 has a low luminosity for a core galaxy, M B = 
-19.6 mag. Michard & Marchal (1994) suspected a peculiar 
asymmetric envelope, but it is difficult to distinguish because 
of nearby stars. It is neighbored by NGC 4319 only 7.4' (56 
kpc at d = 26 Mpc) away, but does not appear to be interacting 
(RC3). Like other core galaxies, it has slightly boxy isophotes 
(100fl4/a = -0.3, BSG), and an excess of X-ray emission above 
that expected from binary stars alone (Pellegrini 1999). 

Our STIS light profile for NGC 4291 (Figure 6) shows a 
slight excess at r < 0."1. No strong evidence exists for a central 
excess in the light profile of the WFPC2 images. 

Our major axis kinematics for NGC 4291 can be compared 
to those of BSG and Jedrzejewski & Schechter (1989). All 
datasets show a centrally rising velocity dispersion which peaks 
at 270-300 km s" 1 . The rotation peaks at -80 km s" 1 at about 5" 
and then declines to about 40 km s" 1 by 12". This overshoot is 
reminiscent of disky galaxies like N3384, but there is no photo- 
metric evidence for a disk in NGC 429 1 . The other authors also 
show minor axis spectra, whereas we have a spectrum offset 
by 29°. NGC 4291 shows very little minor axis rotation (e.g., 
-3.3± 1.2 km s _1 ; Jedrzejewski & Schechter 1989; BSG). The 
velocity profile shows an asymmetric "overshoot" at —5" sug- 
gestive of a co-rotating core. We find lower hA values in our 
ground-based data than do BSG. The BH mass estimate from 
the 2-integral modeling of Magorrian et al. (1998) is 1.9x 10 9 
Mq, while 3-integral modeling applied to the data herein give 
3. 1 x 10 8 M Q (Gebhardt et al. 2003). This is the largest such dis- 
crepancy in our sample. We believe it results from the Magor- 
rian et al. (1998) assumption of isotropy. Wrobel & Herrnstein 
(2000) observed NGC 4291 in 8.5 GHz and placed an upper 
limit on its ADAF accretion rate using the Magorrian et al. 
(1998) BH mass estimate. 

5.6. NGC 4473 

NGC 4473 is an E5 galaxy in the Virgo Cluster. As seen 
in other E5 galaxies, the isophotes are primarily disky (BDM; 



Michard & Marchal 1994). Van den Bosch et al. (1994) note 
some irregularities in the higher-order deviations from ellipses 
between 2" and 4". But ground-based surface photometry gen- 
erally reveals an ordinary, disky elliptical. 

Some searches for an ISM, dust, and emission-line gas have 
found none (Roberts et al. 1991; Michard & Marchal 1994; 
Ho et al. 1997), while Macchetto et al. (1996) detect a small 
amount (1200 Mq) of excited gas in Ha, and Ferrari et al. 
(1999) report a similarly sized dust disk at their limit of de- 
tection. 

Our kinematic measurements are very good for NGC 4473 
because of its relatively high surface brightness (see Fig. 6). In 
particular, the signal in the higher-order moments, h3 and hA, 
is relatively strong (Figure 14). Our ground-based values of h3 
are in good agreement with BSG, but our hA and a values are 
systematically lower. We find hA to be negative at most radii 
and gradually increasing toward the center. hA then shows a 
central dip at STIS resolution. Our a profile shows a gradual 
decrease toward the center where our ground-based values are 
about 10 km s" 1 lower than those of BSG. However, our central, 
Modspec dispersion (183 km s" 1 ) agrees with the average of 1 1 
sources in Hypercat, 179 km s" 1 . Young et al. (1978) published 
long-slit spectroscopy out to 45". They measure a v max of 60 km 
s" 1 and a = 180 km s -1 . Morton & Chevalier (1973) measured a 
higher rotation velocity, lOOkm/s at r= 10". Our ground-based 
Vmax is 75 km s" 1 . Michard & Marchal (1994) point out that this 
rotation is surprisingly small for a disky elliptical. 



NGC 4473 : 




r (arcsec) 



FIG. 15. — Top: deconvolved, WFPC2, V and / surface brightness profiles 
of NGC 4473. Bottom: (V-l) color profile. 

The high-resolution HST imaging reveals NGC 4473 to be 
an unusual galaxy (see Byun et al. 1996 for WFPC1 results, Ta- 
ble 5 for WFPC2 V and / results). Its surface brightness profile 
(Fig. 15) has an asymptotic inner slope, 7 = 0.01 (WFPC2 I- 
band, 7 = 0.014 in WFPC1). This makes it a core galaxy, by 
definition (7 < 0.3; Faber et al. 1997). However, the transition 
between the inner and outer slopes is very gradual (a = 0.66 
for WFPC1; 0.70 for WFPC2), while all the other cores have 
a = 1.0-8.0 (the next smallest being NGC 3608, with the 
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FIG. 16. — Plot of the Nuker-\aw parameters (3 vs. 7 (top), and a vs 7 (bottom) for 72 galaxies. The open circles are the WFPC1 data tabulated in Byun et al. 
(1996). The squares are mostly / band WFPC2 data (V or R if / was not available) from Lauer (2002, private communication). The X 's and filled circles represent 
the 10 galaxies of this paper. For NGC 4473 we show three points: WFPC2 /, V, and WFPC1 V. The four galaxies with very large a values that fall off of the plot 
are marked by circles with tails. 



counter-rotating core). Furthermore, the outer slope of the pro- 
file (J3 = 2.45 and 2.7 for WFPC1 and WFPC2 /-band, respec- 
tively) is very steep - it has the highest j3 among the cores in 
Byun et al. (1996). Figure 16 demonstrates how NGC 4473 
has unusual Nuker parameters. Finally, the absolute magni- 
tude, Mb ~ -19.9 mag, is consistent with power-law galaxies, 
falling faintward of the region of overlap between power laws 
and cores (Faber et al. 1997). Why does a relatively faint, E5, 
disky elliptical have the central profile of a core galaxy? 

A recent merger can provide several explanations for these 
peculiarities. First, the merger could contribute a secondary nu- 
cleus which could artificially flatten the inner surface brightness 
profile and account for the gradual transition between the inner 
and outer slopes. Second, the merger could deposit dust which 
masks the central surface brightness peak. Third, the merger 
may be accompanied by an inflow of gas which forms into a 



stellar disk or torus which, in turn, makes the original peak less 
distinct. 

There is some evidence for each of these scenarios. For ex- 
ample, a disky stellar structure at r < 2.0" is suggested by the 
increasing ellipticity of isophotes towards the center; WFPC1 
surface photometry shows isophotes increasing in ellipticity 
from 0.36 at 15" to 0.5 at 0"5 (van den Bosch et al. 1994). 
Evidence for the recent formation of a stellar disk includes the 
blue color of the stellar disk detected in ground-based images 
(Goudfrooij et al. 1994). The presence of dust is supported 
by the aforementioned observation of central dust and gas (Fer- 
rari et al. 1999). Finally, in support of double nuclei or tori, 
Byun et al. (1996) noted "double nucleus?" for NGC 4473 in 
their Table 1. Perhaps it is no coincidence that NGC 4473 also 
shows the most conspicuous asymmetry in its STIS light profile 
for our entire sample - a shoulder located «0.2" from the peak 



20 




FIG. 17. — WFPC2 greyscale images of the inner part (r <0"7) of NGC 4473. Dark greylevels represent large pixel values in all panels. Top: /-band (left) 
and V-band (right) residuals after subtraction of a model with purely elliptical isophotes. Some isophotes from the ellipse fitting are shown. The center was a 
free parameter. Bottom: /-band (left) and V-band (right) images before model subtraction, displayed with the same orientation and scale. Each has 40 iterations 
of Lucy-Richardson deconvolution. Minimally-smoothed contours are overlayed with linearly-spaced levels. The inner isophotes become so boxy as to appear 
"peanuf'-shaped. 



(see §4.5, Fig. 6). 

The above conclusions that dust, blue colors, and secondary 
nuclei are present are all based on only marginal detections. 
Thus, we have inspected the deconvolved HST WFPC2 images 
for corroborating evidence. The / and V band images (Figure 
17) have 2000 s and 1800 s total exposure, respectively. They 
provide higher resolution than the WFPC1 data used by van den 
Bosch et al. (1994). Our HST (V-I) color profile in Figure 15 
does not confirm the blue color of the disk, although blue and 
red light may be mixed within each elliptical isophote. A dou- 
ble nucleus is not present in the form of two, distinct maxima. 
However, the contours appear compressed along the minor axis 
so as to make isophotes "peanut" shaped at r 0.2". Ellipse 
fitting reveals that the strong diskiness in isophotes at r > 1 .6" 
(peaking at 100a4/a = 2.3 ±0.5 atr=4.3") suddenly gives way 
to very strong boxiness at r < 1.6"; the lOOa^/a drops below 
-5.0 at r =0.23" and then returns to values consistent with at 
the resolution limit (r « 0. 1"). The noisiness of the deconvolu- 
tion process should be taken into consideration; the peak pixel 
is not in the same position in the two bottom images and the 
minimally-smoothed contours are noisy. However, the position 
of the boxiest isophotes is robust, occurring at r w 0.2" in both 
the V and / band deconvolved images and reaching magnitudes 
as low as 100a4/a = -5 in both datasets. The top panels of Fig- 
ure 17 shows that subtracting a purely elliptical model leaves 
behind an " x "-shaped residual in both cases. The offset of the 
peak residual from the best-fit ellipse centroids is also robust, 
suggesting an asymmetric enhancement about 0.1" from those 



centroids. At r > 1" (not shown), the residual pattern inverts 
to a dark "+", indicating the diskiness previously reported by 
others. 




FIG. 18. — V-band WFPC2 data shown with the STIS slit overlayed to match 
the observation. The circle is 0"2 wide. North is 173° clockwise from up, and 
the position along the slit used elsewhere in this paper increases towards the 
left. 

The asymmetry we see in the residuals after subtracting an 
elliptical model are one way to explain the shoulder in the light 
profile along the STIS slit. A more direct test is to simulate the 



21 



STIS slit aperture on the WFPC2 data and create light profiles. 
These profiles are shown in Fig. 6. The light is more concen- 
trated near the center because of the greater spatial resolution in 
the WFPC2 data, especially in the deconvolved data. Neverthe- 
less, the profiles are higher to the right of the center than to the 
left under the "double nucleus" label. Figure 18 is an overlay of 
the 0."2 STIS slit onto the V-band WFPC2 image. Notice that 
the STIS slit is skewed with respect to the major axis by 20°, 
more than any other galaxy in our sample (this was necessary 
to acquire the guide stars). Nevertheless, the position of the slit 
allows the passage of some light from the region which created 
the " x "-shaped residuals. 

There are no obvious signs of a recent merger in the kinemat- 
ics. No counter-rotation is seen in the velocity profile, and the 
velocity "overshoot" seen at 5" is typical of ellipticals with em- 
bedded disks. The photometry and STIS light profile (Fig. 6) 
suggest a possible secondary nucleus at r « 0.2", so a kinematic 
deviation might be expected at STIS resolution. The unsym- 
metrized STIS dispersion in Figs. 14 and 7 does, in fact, show 
a smooth gradient across the center, but it is only marginally 
significant. NGC 4473 , like NGC 3384, has a lower central 
dispersion in the STIS data than in the ground-based data, in- 
dicating a kinematically cold stellar population near the center. 
These two galaxies are interesting to contrast: both show kine- 
matic and photometric evidence for a cold stellar disk all the 
way into the center, yet NGC 4473 has much lower (v/a) than 
NGC 3384. Moreover, NGC 4473 has a core profile with an 
indistinct center, while NGC 3384 has a power-law profile with 
the strongest central peak in our sample. 

We conclude that NGC 4473 is a peculiarity - a disky galaxy 
lacking a strong central cusp - and that its peculiarities prob- 
ably originated with a merger. The occurrence rate of such 
'disky cores' is low, judging by the study of Rest et al. (2001) 
which classifies only 1 of 9 cores as "disky" (within a sample 
of 57 galaxies). Any obvious evidence for a recent merger is 
absent; there is no counter-rotating stellar or gas system, no 
erratic dust, and no multiple nuclei. But current work on the 
formation of cores (e.g., Faber et al. . 1997; Milosavljevic & 
Merritt 2001; Makino 1997) involves the merger and coales- 
cence of two galactic nuclei with BH, and these scenarios pre- 
dict that observability of the merger will diminish continuously. 
Some cores have been identified with easily resolved double 
nuclei (e.g., NGC 4486B, Lauer et al. 1996), and recently, six 
cores have been identified with subtle "central inversions" in 
their surface brightness profile which are likely to be toroidal 
stellar systems (Lauer et al. 2002). NGC 4473 may be another 
case of subtle central structures resulting from coelescing bi- 
nary nuclei. Like most of those in Lauer et al. (2002), it has 
a poorly defined central peak and "peanut"-shaped isophotes. 
NGC 4473 differs in that it has disky isophotes at r >l". Alter- 
natively, no BHs were contributed by the last merger, only gas 
which has subsequently formed a substantial disk at r > 1 .5"and 
a non-cuspy structure (possibly a torus) near the limits of reso- 
lution. Future imaging at slightly higher resolution should clar- 
ify the central structure. 

5.7. NGC 4564 

NGC 4564 is classified by Sandage & Tammann (1980) as 
an E6, but it has characteristics of an SO. First, its isophotes 
show disky deviations from an ellipse at large radii. The pa- 
rameter an/a begins to climb at 10", reaching about 0.025 at 
20" and then tapering down to 0.01 (BDM). Such diskiness at 



intermediate radii is also seen in NGC 821. Second, its elliptic - 
ity increases gradually from 0.2 at r = 3" to 0.6 at 30" (BDM). 
Sahu et al. (1996) decompose their R band image into an el- 
liptical bulge (obeying a de Vaucouleurs law) and a disk. Their 
scale length for the disk is r s = 28.14" with axial ratio b/a = 0.3. 
Third, Sil'chenko (1997) finds a younger stellar population in 
the center of the disk (8 Gyr) compared to the rest of the galaxy. 

High-resolution images were obtained with WFPC1 (PID 
2607, Byun et al. 1996) and WFPC2 (Jaffe, PID 6357; Faber 
et al. 1997). In the center, NGC 4564 appears more bulge- 
dominated, with rounder isophotes and boxy deviations from 
ellipses. For example, van den Bosch et al. (1994) find that 
the ellipticity is only 0.16 at 0"5 and that a^ja is significantly 
negative around 2"0. They also find sharp variations in PA and 
e near 2."0. This is indicative of dust; however, Roberts et al. 
(1991) do not find evidence for an ISM. 

Our ground-based kinematics for NGC 4564 (Figure 19) 
can be compared to those of BSG and Halliday et al. (2001). 
Our spectra have relatively low S/N (Figure 3). The pub- 
lished profiles are in good agreement with ours, but are plot- 
ted to greater radii (r = 35"and 40", respectively). The ve- 
locity dispersion declines by a factor of 2 from 0"to 10", so 
one should expect large variations in published values of a- 
dependent quantities. For example, Bender 1988 give a v max /a 
of 1.11 ( Vmax — 150±5 km s a— 135 ±5 km s 1 ) while we have 
1.3 ±0.2 0w=147±4kms- 1 , a= 112±17kms _1 , see Table 
8 ). The average of 6 central velocity dispersions in Hyper- 
cat is 158 km s" 1 , while we measure ctq ~170 km s -1 in all of 
our Modspec PAs. Asymmetric (h3) deviations of the LOS VDs 
from a Gaussian are not significant, but hA appears to be nega- 
tive at r < 10". Halliday et al. (2001) show h.3 becoming nega- 
tive beyond 10", coincident with an increase of 04 above 0. The 
sudden changes in the PA, ellipticity, cos 36*, and cos4# coeffi- 
cients at 2"0 (van den Bosch 1994) are coincident with subtle 
bends in our profiles of symmetrized v, a, and hA. 

The STIS kinematics feature a central velocity dispersion 
peak of 220±20 km s" 1 , well above the ground-based result. 
The higher resolution also reveals rotation of ^50 km s" 1 in 
the inner 1 "0. h3 and hA are fairly noisy, but the symmetrized 
curves indicate a gradient in h3 across the inner 0"5, and pre- 
dominantly negative hA. The kinematics allow a secure detec- 
tion of a BH (Gebhardt et al. 2003), but little evidence for nu- 
clear activity is found in X-rays (Pellegrini 1999) or radio con- 
tinuum (Wrobel & Herrnstein 2000). 

5.8. NGC 4649 

NGC 4649 (M60) is a giant elliptical in the Virgo cluster, 
comparable in luminosity to M87. NGC 4649 stands out in our 
sample as the galaxy with the lowest central surface brightness 
(Fig. 6) and highest luminosity. It is the only galaxy with a 
large enough dispersion so that the crowding of the Ca II triplet 
absorption lines becomes problematic. This will be discussed 
further in §6. 

We find a 10 9 3 M BH in NGC 4649 (Gebhardt et al. 2003), 
only about 30% less massive than the BH in M87 (Macchetto 
et al. 1997). Fortunately, it has far weaker nuclear activity than 
M87, thus simplifying a BH search. Nuclear activity is detected 
in the radio (Fabbiano et al. 1987) and but not in the optical or 
X-rays (Byun et al. 1996; Di Matteo & Fabian 1997). Stanger & 
Warwick (1986) also find only low-level extended radio emis- 
sion out to ~3."0. Di Matteo & Fabian (1997) estimate that the 
total measured core flux from radio to X-ray is lower than 10 41 
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FIG. 19. — Kinematic profiles for NGC 4564 (see Fig. 9 for the meaning of the symbols). For NGC 4564, all ground-based PAs are from Ca spectra. 



erg S . 

NGC 4649 is separated by only 2.5' (~ 12 kpc for a distance 
of 16 Mpc) from the spiral NGC 4647. Many studies report 
that they are non-interacting (e.g., Sandage & Bedke 1994), 
although the spiral is clearly asymmetric (see Koopman et al. 
2001). It is given an E2 morphology in the RC2 (de Vau- 
couleurs et al. 1976), but Sandage & Bedke (1994) think it is an 
S0i because of its "prominent extended envelope." NGC 4649 
has an extended X-ray halo (Bohringer et al. 2000), and it has 
been identified with its own compact subgroup within the Virgo 
cluster (Mamon 1989). 

CCD surface photometry has been performed by Peletier 
et al. (1990), Caon et al. (1990), and Michard & Marchal 
(1994). There are no significant deviations from ellipses in 
the isophotes between l"and 8", and then they become mildly 
boxy (100fl4/a w -0.5). The maximum ellipticity measured by 
Peletier et al. (1990) is « 0.2 which suggests a nearly face- 
on inclination (i < 36 °). This makes difficult the detection of 



disky isophotes (Michard & Marchal 1994). Instead, support 
for the SO classification is found in the surface brightness pro- 
files which show an outer envelope above the r 1,/4 -law bulge. 
There is no strong evidence for a tidal influence from nearby 
NGC 4647. HST WFPC2 V-band photometry is discussed by 
Byun et al. (1996) and Faber et al. (1997). NGC 4649 has a 
core profile with a break radius of 3.58". 

Our ground-based kinematics can be compared to Fisher 
et al. (1995; hereafter FIF), De Bruyne et al. (2001), and BSG 
FIF report an average rotation of 87 km s" 1 , and their rotation 
curve climbs from 60 to 120 km s" 1 between r =20 to 60". Sim- 
ilarly, our rotation curve is climbing at 45" where it is ~110 
km s" 1 (symmetrized). In general, NGC 4649 has strong ro- 
tational support compared to other giant ellipticals. Its large 
(v/cr)* stands apart from the brightest cluster galaxies in the 
study of FIF. Our unsymmetrized rotation curve is also fairly 
asymmetric with respect to the galaxy center. The data of De 
Bruyne et al. (2001) demonstrate this even better with a rotation 
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FIG. 20. — Kinematic profiles for NGC 4649 (see Fig. 9 for the meaning of the symbols). For NGC 4649, all of the ground-based data are from Mg spectra except 
the dashed line (PA = 22°) which is Ca spectra. 



curve out to 90", showing an ~ 70 km s" 1 maximum difference 
in rotation between the two sides. This asymmetry along the 
major axis provides some support for the idea of an interaction 
with NGC 4647. 

Our velocity dispersion profile is consistent with that of BSG 
and De Bruyne et al. , but FIF has systematically lower values. 
Most dispersion profiles peak at ~ 400 km s" 1 and then rapidly 
decline to ~ 300 km s" 1 by r = 20"; they show nearly a factor of 
two decrease from the center to the effective radius (r e « 80"). 
De Bruyne measures positive hA values within 10" whereas we 
measure negative /i4 here at all of our PAs. 

The STIS Ca II triplet data reveal an enormous central veloc- 
ity dispersion. The entire inner 1 "0 is above 450 km s" 1 , and it 
appears to rise to over 600 km s" 1 . However, the uncertainties 
are high because of the ill-defined absorption lines (§6). This 
differs from the case of IC 1459 (Cappellari et al. 2002), an- 
other giant elliptical (Mb ~ -21 .4), which has ctq 350 km s -1 



in both STIS and ground-based observations. Another unusual 
(but uncertain) feature is the strong rotation within 1 ."0. The 
maximum STIS rotation, ^150 km s _1 , exceeds the maximum 
ground-based rotation. 

5.9. NGC 4697 

This is an E6 galaxy with a stellar disk along the apparent 
major axis (Carter 1987; Goudfrooij et al. 1994a ). NGC 4697 
has more dust and gas than the other galaxies in this sample. 
Molecular gas was detected by Sofue & Wakamatsu (1993), and 
the high IRAS 100 /Ltm flux densities also imply cool gas. The 
radial color gradient is larger than the metallicity gradient from 
the Mg2 index (Peletier 1990). HST WFPC1 imaging reveals an 
organized dust disk (Lauer et al. 1995). This dust was not suc- 
cessfully imaged from the ground (Veron-Cetty & Veron 1988; 
Kim 1989), but was suggested by the red nucleus in B-I (Goud- 
frooij et al. 1994a). The Ha+[N II] emission proved difficult 
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FIG. 21. — Kinematic profiles for NGC 4697 (see Fig. 9 for the meaning of the symbols). For NGC 4697, all ground-based PAs are from Ca spectra. In the top 
left subpanel, we include filled circles for [N II] emission line kinematics (see Pinkney et al. (2003)). 



to detect (Trinchieri & di Serego Alighieri 1991; Kim 1989). 
Long-slit spectroscopy by Goudfrooij et al. (1994b) gives an 
Ha+[N II] flux of 2.4xl0" 13 erg s" 1 cm" 2 . The narrowband 
imaging by these authors shows Ha+[N II] emission extending 
out to 35". Ground-based imaging (e.g., Peletier 1990) indi- 
cates embedded, disky isophotes. The parameter lOOa^/a is 
over 1.8 until r= 20"and then it declines to 0. 

A great deal of spectroscopy exists for comparison with our 
Modspec results. Bertola & Capaccioli (1975) first pointed out 
that ellipticals are not rotationally supported using NGC 4697. 
It has also been used for case studies showing that 3-integral 
model are an improvement over 2-integral models (Dejonghe et 
al 1996; Binney et al. 1990). 

Carter (1987) points out that there is a large variation in v max 
measurements for this galaxy. Bertola & Capaccioli (1975) 
find 65 km s , Illingworth (1977) gets 90 km s -1 , and Davies 
(1981) derives 146 km s" 1 . We measure a ground-based v maK = 



101 ± 13, and a STIS v max of 110 ±21 km s _1 . A likely ex- 
planation is the presence of a cold, rapidly rotating, stellar disk 
(composing more than 10% of the light along the major axis) 
which has a varying effect on the measured radial velocity for 
different methods. Figure 21 includes the rotation curve of the 
excited gas (black points) measured from [N II] emission with 
STIS (Pinkney et al. 2003). This provides an interesting con- 
trast to the stars in that the gas shows less pressure support be- 
yond 0."2. 

The mean central velocity dispersion given by Hypercat is 
174 km s" 1 . This is a bit higher than our central, GB, velocity 
dispersions of 157, 159, and 164 km s" 1 , for our 0, 48, and 70° 
slit orientations, respectively. We measure a higher central dis- 
persion with STIS, 171 km s -1 . This galaxy stands out as one 
of the few with a local minimum at r = in the STIS velocity 
dispersion profile (Fig. 21). It is also unusual in that its h4 
parameter is significantly positive (see Fig. 25) 
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FIG. 22. — Kinematic profiles for NGC 7457 (see Fig. 9 for the meaning of the symbols). For NGC 7457, all ground-based PAs are from Ca spectra. 



5.10. NGC 7457 

NGC 7457 is the smallest, lowest luminosity galaxy in our 
sample. It is an SO galaxy, like N3384. NGC 7457 is another 
example of a pseudobulge (see §5.3) in our sample, based on 
photometry and kinematics. Tomita et al. (2000) find a point- 
like, bluish nucleus in the WFPC2 (V,I) data, and Peletier et al. 
(1999) find a young stellar population in the bulge. Michard & 
Marchal (1994) detect a small bar. 

Lauer et al. (1991) report on the WFPC surface brightness 
profile and find it unresolved and steep. The improved HST 
WFPC2 imaging reveals a surface brightness profile with a very 
shallow outer power law; its j3 = 1 .05 (Lauer 2002, private com- 
munication) is the smallest in a sample of 29 early types. How- 
ever, the surface brightness profile in Andredakis et al. (1995) 
has a Sersic parameter n w 6.2 which is higher than the values 
n <2.0 more typical of pseudobulges. Most impressive is the 
strong surface brightness peak in deconvolved WFPC2 F555W 



data (Lauer 2002, private communication) which jumps 4 mag- 
nitudes in the inner Of'l. Tomita et al. (2000) do not interpret 
this as an AGN because of the lack of Ha and H/3 emission re- 
ported by Ho et al. (1997). However, Gebhardt et al. (2003) find 
a central drop in calcium triplet equivalent widths and interpret 
the peak as AGN light. 

We find evidence for pseudobulge-like kinematics for NGC 
7457. Our Figure 23 demonstrates that the velocity dispersion 
falls low on the Faber-Jackson relation. This result is robust for 
different global values of a from Table 8. Our v max la values 
(Table 8) are not especially large, however, given the ellipticity. 

6. DISCUSSION 

6.1. Central Velocity Dispersions 

To what extent does the STIS central velocity dispersion 
drive the detection of the nuclear black holes? Modeling is re- 
quired to fully answer this question; however, a strong peak 
in dispersion in the central 0"5 (~ 45 pc at d = 17 Mpc) is 
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difficult to explain without a central dark mass (van der Marel 
1994b). In Figure 23, we compare the central STIS dispersion 
to our ground-based, "effective" dispersion, u e (i.e., the rms 
velocity relative to the systemic velocity of the galaxy, aver- 
aged over a slit width of 1" extending to R e \ GOO). In 7 out 
of 10 cases, the STIS central value is higher. The most ex- 
treme differences are NGC 4649 (162 km s" 1 , 42%), NGC 3608 
(114 km s"\ 55%), NGC 4564 (58 km s" 1 , 37%), NGC 4291 
(64 km s~\ 17%), and NGC 7457 (38 km s" 1 , 60 %). Three 
of these have core profiles, and it appears that the core galaxies 
show the largest differences. NGC 7457 has a 60% increase 
from u e to a S Tis, the greatest fractional rise in the entire sam- 
ple. Figure 23b demonstrates that the fractional increase does 
not correlate with increasing dispersion. Likewise, there is no 
correlation with M B Mi ge (not shown). NGC 4473, unlike the 



other core galaxies, has a ustis lower than a e . This adds to the 
growing list of differences between NGC 4473 and the cores 
(§5.6). The other galaxies with lower vstis are NGC 3384 and 
4697, which are also disky. However, NGC 4697 has a higher 
(Jstis than our alternative ground-based dispersions oo.gb and a 
in Table 8. Eight of our galaxies, including NGC 3384, show a 
local maximum at r =0"0 in their STIS dispersion profile. NGC 
4473 shows only a slight central, local maximum, while NGC 
4649 and NGC 4697 show local minima in their symmetrized 
a profiles. Thus, the majority of galaxies (70 - 80%) show a 
central peak at STIS resolution, and about 70% have a larger 
o"o,sr/s than a e , Only one galaxy, NGC 4697, does not obey 
either criterion. 

In Figure 23c, we show the Faber-Jackson (1976) relation for 
our 10 galaxies. The a e is an adequate surrogate for the tradi- 
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tional central dispersions measured through ~ 2 x 2" apertures, 
but use of (Jstis adds scatter to the relation presumably because 
of the strong dispersion gradient induced by the central black 
hole. We also find a marginal correlation between the ratio 
o'STis/< J e and the ratio v ma */<7, where v max is the maximum ro- 
tation in the ground-based data and 7r is the average dispersion 
in Table 8. Apparently, galaxies with more rotational support 
have less of a dispersion peak at STIS resolution. The proba- 
bility that no correlation exists is < 3% (using Cox regression 
and Kendall's Tau, Isobe et al. 1986) This correlation is present 
but not as strong if v max /a is replaced by (v max /a)* . 

6.2. Low Surface Brightness Giant: NGC 4649 

The galaxies in the BH literature, including our own sample, 
show a preference for high surface brightness in core galaxies 
since these are easier to observe. In fact, the measurement of 
dispersion profiles in core galaxies is so difficult that most of 
the BH detections in these galaxies rely on emission lines from 
rotating gas disks [e.g., M87 (Harms 1994), NGC 4261 (Fer- 
rarese et al. 1996), M84 (Bower et al. 1997)]. In the case of 
the giant E3, IC1459, Cappellari et al. (2002) have measured 



the BH mass using stellar kinematics as well as gas kinemat- 
ics. However, this galaxy does not have an especially faint sur- 
face brightness, nvfi.i" — 15.3 (Carollo et al. 1997). Ground- 
based stellar kinematics have provided a BH mass estimate for 
M87 (Dressier & Richstone 1990; van der Marel 1994a) but the 
BH was not required to fit these data. The preference for high 
surface brightness cores could conceivably introduce a bias in 
plots of Mbh -cor Mbh - Lb . Moreover, the measurement of BH 
mass using gas kinematics gives a different answer than stellar 
kinematics in the only giant where they have been compared 
(IC1459, Cappellari et al. 2002). Thus, NGC 4649 is important 
as the only BH measurement from stellar kinematics in a galaxy 
which is representative of low-surface brightness cores. 

Not surprisingly, we find that the LOSVD is particularly dif- 
ficult to measure in NGC 4649. We requested sufficient expo- 
sure time to measure Ca lines in a galaxy with n v(j n x =\5.9 mag 

arcsec" 2 and a « 350 km s" 1 . However, this galaxy appears to 
have a very large central velocity dispersion at STIS resolution. 
This broadens the lines to such an extent that their detection is 
difficult even at this S/N (Figure 24). Our MPL technique was 
able to measure an LOSVD with a >500 km s -1 . Naturally, the 
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upper limits are more weakly constrained than the lower limits. 
The FCQ technique gave more erratic results (not shown in Fig. 
20) for these data. 

6.3. Gauss-Hermite Parameters 

Our parameterized stellar LOSVDs suggest properties of a 
galaxy's orbital structure that are elucidated by the 3-integral 
modelling. This merits a comparison between the two. Geb- 
hardt et al. (2003) find the tangential component of the velocity 
dispersion tensor to be stronger in the central bins of the galaxy 
than in the bins near R e /A (see Figure 25b). The hA parameter 
also indicates a tangential bias when it takes a negative value 
(i.e., the LOSVD is boxy). Indeed, we find that our ground- 
based hA profiles usually have a central dip (NGC 821 is a good 
example). Moreover, the central, STIS-measured hA parame- 
ters, although noisier, are also predominantly negative. To mea- 
sure a more reliable central STIS hA, we averaged the hA values 
from the innermost 3 bins. Figure 25 demonstrates that these 
values of hA are negative in 9 out of 10 cases. Individual cases 
are not significantly negative, but altogether they suggest a real 
trend. NGC 4473 and NGC 3608 are the most significantly 



negative. There is one outlier, NGC 4697, which has a positive 
hA and this pulls the weighted mean to -0.027 ±0.01. When 
this point is excluded, the weighted mean is -0.05 ± 0.01. We 
looked for correlations between the central hA value and other 
photometric and kinematic parameters. The large uncertainties 
in the hA values and the small sample do not permit any signif- 
icant results. A positive correlation is hinted at by Figures 25a, 
and 25c. 

It is puzzling that the disky elliptical, NGC 4697, shows posi- 
tive central hA values while the other disky galaxies (NGC 821, 
3384, 4564, and 4473) have negative central hA values. One 
might search for other ways in which NGC 4697 differs from 
the other four. The most obvious difference is its prominant 
dust disk at r < 3.5". This disk is inclined at 77° and has a 
sharp outer cut-off and a less clear-cut inner radius (Pinkney 
et al. 2003). One would expect the same LOSVD regardless 
of the presence of the dust disk, if the dust disk is an infinitely 
thin screen centered within an axisymmetric galaxy. However, 
one could invoke a disk with finite thickness to preferentially 
obscure, say, an equatorial cold disk component and thereby in- 
fluence the LOSVDs. A second possible scapegoat is the stellar 
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disk itself, which may be more (or less) prominent in the cen- 
tral STIS bins of NGC 4697 than the other four disky galaxies. 
One expects a superposition of two LOSVDs with the same 
centroid but different widths to produce an LOSVD with pos- 
itive h4. Thus, we may find positive values of h4 in galaxies 
where a bulge and disk population overlap. However, there is 
little evidence in ground-based kinematics that NGC 4697 has 
a significantly different disk contribution than the other galax- 
ies. In fact, NGC 4473 seems to have the strongest contribution 
of light from a stellar disk at r < 5" : it was the only galaxy 
that warranted the inclusion of a disk component in the mod- 
eling (Gebhardt et al. 2003). Finally, we see in Figure 25 that 
NGC 4697 has the lowest mean h3 at r > 2" of the entire sam- 
ple. This is another artifact of diskiness, or more correctly, of 
strong rotation. Unfortunately, the kinematic properties mea- 
sured at r > 1 "0 do not necessarily predict the properties at 
r < 0"1 where the h4 values in question are measured. 

The bottom three panels of Figure 25 complement the inves- 
tigation of Gauss-Hermite parameters. First, Figure 25d shows 
our 10 galaxies plotted on a correlation found by BSG Our 
galaxies do not follow the trend well, but this is mostly because 
of NGC 4564. The v max la in this figure uses a, but if we used 
cq.gb or a e (Table 8), v max /a would be 0.9 for NGC 4564 instead 
of 1.3. NGC 4564 happens to have the most extreme variation 
between a and our other two ground-based a estimates for rea- 
sons that are clear in Figure 19. Second, Figure 25e demon- 
strates that galaxies with strong rotational support tend to be 
disky. This has been established with larger samples by BSG. 
Third, Figure 25f shows a correlation between the asymptotic 
inner slope 7 of the best-fit Nuker-law surface brightness pro- 
file, and the deviations of isophotes from ellipses. This follows 
from the results of, e.g., Faber et al. (1997) wherein the disky 
galaxies (1 0004/0 > 0) tend to have power-law profiles and the 
boxy-isophote galaxies (lOOa/t/a < 0) tend to have core pro- 
files. This figure is intended to underscore the peculiarity of 
NGC 4473 as a "disky core. " 

6.4. Demographic Results 

Our sample of early-type galaxies includes a wide range of 
velocity dispersions (70 to 385 km s" 1 ) and BH masses (10 6,5 to 
10 9 3 Mq), allowing us to address questions of BH demograph- 
ics (§1). A word should be said about selection effects near 
the high-mass end of the M/,h - <J relationship. The measure- 
ment of LOSVDs becomes increasingly more time consuming 
with decreasing surface brightness. Therefore, a selection bias 
is introduced whereby, given two core galaxies of the same lu- 
minosity, or the same dispersion, one is compelled to select the 
one with the higher surface brightness. We selected one chal- 
lenging galaxy located in this region of parameter space, NGC 
4649, with /no,v = 15.9 mag arcsec" 2 . Nevertheless, none of 
our cores, including NGC 4649, represent the lowest surface 
brightnesses known in their luminosity range. It is not certain 
whether the lowest surface brightness galaxies obey the same 
demographic trends. 

With the discovery of the Mbh - o~ relation, it is natural to look 
for exceptions to the rule. It appears that the relation holds for 
Sa and Sb galaxies with classical bulges, and even later types 
(Sc - Sd) with pseudobulges (Kormendy et al. 2002). The Sc 
galaxy M33 appears exceptional in that it has a BH mass upper 
limit significantly lower than that predicted using the dispersion 
of its nuclear star cluster (Gebhardt et al. 2001). This galaxy, 
however, lacks a hot spheroidal component and so it is debat- 



able whether the Mbh — & relation is applicable. 

In Figure 26, our ten BH masses are plotted against velocity 
dispersion. For just these galaxies, we fit an - a relation of 
the form 

log lo (M M /M ) = a + /31og 10 (cr e /200 km s" 1 ). (3) 

The fitting procedure is identical to the one advocated by 
Tremaine et al. (2002), where we have assumed 5% error in 
<r e , and the uncertainty in Mbh is the measured one combined 
with an extra 0.34 dex of intrinsic dispersion, which gives a x 2 
per degree of freedom of unity. We find a = 8.06 ±0.13 and 
j3 = 3.67 ± 0.70. Note that the BH masses for the galaxies used 
in this fit have been refined since GOO. All galaxy distances are 
now taken from Tonry et al. (2001), and more dynamical mod- 
els have been run to better probe chi-square space (Gebhardt 
et al. 2003). Nevertheless, the best fit is nearly indistinguish- 
able from the GOO fit. All of these same refinements to our 
sample are incorporated into the preferred fit by Tremaine et al. 
(2002), which includes an additional 21 galaxies from the lit- 
erature. The most deviant point is for NGC 2778. This faint 
galaxy has a lower S/N than most. Also, its sphere of influence 
is unresolved. Consequently, NGC 2778 has the least confident 
BH detection when modeled (Gebhardt et al. 2003). 
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FIG. 26. — Updated plot of M^-cr for only the galaxies in this paper (la- 
belled squares). The best-fit line to these 10 galaxies is shown as a solid line. 
We overlay the original fit by Gebhardt et al. 2000 (dotted line), and the pre- 
ferred fit by Tremaine et al. 2002 (dashed line). 

In hierarchical scenarios, the formation of bulges involves 
mergers which are believed to have occurred during the first 
few Gyr of the universe. But what is the effect of a merger oc- 
curring < 1 Gyr before the present? Will the likely increase in 
bulge dispersion be accompanied by the appropriate increase in 
BH mass? We have inspected our galaxies for merger candi- 
dates to see if they are systematically displaced on the Mbh - & 
relation. Three of our galaxies show marginal evidence for re- 
cent mergers: NGC 3608, 4473, and 4697. NGC 3608 has a 
counter-rotating core and patchy dust. NGC 4473 has unusual 
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photometric properties that may be related to a merging com- 
ponent (§5.6). NGC 4697 has a dust ring which may have been 
deposited by a merger. As can be seen in Figure 26, all three 
of these fall reasonably close to the M hh - a relation. Two of 
these, NGC 4697 and 4473, are among our most secure BH 
detections. 

7. CONCLUSIONS 

We have presented long-slit spectroscopy from STIS aboard 
HST and the MDM 2.4-m telescope that is used to derive black- 
hole masses for 10 galaxies by Gebhardt et al. (2000). Ninety 
percent of the galaxies have a centrally peaked velocity dis- 
persion and/or a higher dispersion from STIS than the ground- 
based Modspec. The galaxies with the strongest rotational sup- 
port, as quantified by vmax/&stis, have the smallest dispersion 
excess at STIS resolution. When only our ten galaxies are used 
for a fit to Mbh-cr, the slope is 3.68, which is similar to the value 
of 3.75 found by GOO for 26 galaxies. We review individual 
galaxies and identify candidates for recent merger: NGC 4697, 
3608, and 4473. These galaxies are not outliers in the M^-cr 
relation. We also identify pseudobulges in NGC 3384 and 7457 
and find that they also obey the Mbh - a relation. Finally, we 
find a trend toward flat-topped line-of-sight velocity distribu- 
tions (i.e., negative values of h4) at the center of the galaxies, 
implying a tangential bias in the stellar velocity dispersion. 
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Table 1 
Target Galaxies 



JName 


type" 


M B " 


Frofile^ 


Uist. d 


Scale^ 


Hv 1 


PA s 


i e 


NGC 


mag 


Type 


Mpc 


pc/" 


mag/D" 


deg 


II 


821 


E67 


-20.41 


"A 


24.1 


117 


14.3 


25 


50 


2778 


E 


-18.59 


\ 


22.9 


111 


14.9 


40 


16 


3384' 


SBsO- 


-18.99 


\ 


11.6 


56 


13.7 


53 


25 


3608 


E2 


-19.86 


n 


23.0 


111 


14.7 


75 


34 


4291 


E 


-19.63 


n 


26.2 


127 


14.9 


110 


17 


4473 


E5 


-19.89 


\n 


15.7 


76 


14.7 


100 


26 


4564 


E 


-18.92 


\ 


15.0 


73 


14.1 


47 


20 


4649 


E2 


-21.29 


n 


15.9 


88 


15.9 


105 


69 


4697 


E6 


-20.24 




11.7 


57 


14.4 


70 


72 


7457 


SArsO-? 


-17.69 




13.2 


64 


16.2 


130 


32 



a Galaxy type from NED. 

b Blue absolute magnitude, M B , calculated using surface brightness fluctuation 
(SBF) distances. 

c Surface brightness profile type characterized as either \ (power-law) or n (core). 
Sources are Faber et al. (1997), Byun et al. (1996, for NGC 4473), and Lauer (2002, 
private communication, for NGC 821, 2778, 4291). 

d Distances are from SBF measurements in Tonry et al. (2000). 

e Spatial scale derived from col. 5. 

f Central "fiy" is the V-band surface brightness measured at r =0"1 . 
g Position angle from the RC3. 

h Effective radius derived from the RC3, in arcseconds. 
iNGC 3384 also identified as NGC 3371. 
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Table 2 

Long-Slit Spectrograph Configurations 



Code 


Instrument 
+ grating 


Slit size 

"x" 


— \ :r 

^cen 

A 


A-range " 
A 


Dispersion 
A pix" 1 


1 b (Jomp. line a * 

A ( km s" 1 ) 


Spatial scale J 
" pix" 1 


1 


STIS + G750M 


52x0.1 


8561 


82/5-884/ 


1.1089 


0.45 (17.5) 


0.050/1 


2 


STIS + G750M 


52x0.2 


8561 


8275-8847 


1.1089 


0.68 (23.8) 


0.05071 


3 


STIS + G750M 


52x0.1 


6581 


6295-6867 


1.1130 


0.37(16.9) 


0.05071 


4 


Charlotte+831g/mm 


760x0.8 


8500 


7801-9199 


1.44 


0.98 (35.) 


0.593 


5 


Wilbur+831g/mm 


760x0.9 


8500 


7580-9420 


0.909 


1.10(39.) 


0.371 


6 


Charlotte+ 1 200g/mm 


760x0.9 


5175 


4689-5661 


1.11 


0.75 (37) 


0.593 



"Central wavelength and wavelength range. STIS values taken from Leitherer et al. (2001, pp. 231, 234). 

b Reciprocal dispersion was measured using our own wavelength solutions. The distribution of dispersion solutions 
for a given dataset had a a w 0.00015 Apix -1 . The average dispersion given in the Handbook for G750M is 1 . 12 Apix -1 . 

Tnstrumental line widths measured by fitting Gaussians to emission lines on comparison lamp exposures. This 
gives an estimate of the instrumental line width for extended sources. We use ^5 lines per exposure, and at least 5 
measurements per line. In the cases of the G750M setups, the comparison lamp exposures were unbinned, while the 
galaxy spectra were binned by 2X along the dispersion axis. The 2X binning increases the measured widths by, e.g., 
~25% for the 0."1 slit, and by ~ 3% for the 0."2 slit at 8561 A. Leitherer et al. (2001, p. 300) give these instrumental 
line widths for point sources: er=13.3, 15.0, 16.7 km s" 1 for the first three G750M setups in this table. 

d The spatial scale along the slit is constant, but it varies across the slit from one grating to the next. It is 0.05597 " 
pix" 1 for G750M at 8561 A and 0.05465 "pix" 1 for G750M at 6581 A (STIS ISR 98-23). 

Note. — Most of the STIS setups used a CCD readout with 2X overbinning along the dispersion, and a readnoise of 
~ 1 e"pix _1 . The Modspec setup values are given by the program "modset" by J. Thorstensen. The slit widths, spectral 
and spatial resolutions varied for the Modspec observations; characteristic values are shown here . 



Table 3 
STIS Observations 



IName 


Lxposure 


Ubs. Date - 


PA U 


APA L 


Setup d 


NGC 


seconds 


DDMMYY 


degrees 


degrees 


Code 


821 


27478 


131098,080999 


-150 


5 


1 


2778 


15026 


030699 


50 


10 


1 


3384 


27744 


100598,160499 


58 


5 


1 


3608 


12950 


140598 


94 


19 


2 


4291 


16013 


210498 


110 





2 


4473 


12840 


041297 


-100 


-20 


2 


4564 


15060 


080398 


-143 


-10 


1 


4649 


46467 


090699 


95 


-10 


2 


4697 


24820 


270799,200899 


68 


-2 


1,3 


7457 


15282 


240999 


130 





1 


HR7576 e 


5x3.2s,5x3.2s 


140599 


-173 




1,2 


HR6770 


5x2.2s,5x3.2s 


060999 


54 




1,2 



a The date at the end of the visit. 

b Position angle of the STIS slit. From the FITS headers, 
difference between the STIS position angle and the RC3 PA. 
d Code for the spectrograph setup used (see Table 2). 

template stars were HR7576 (V=5.03 mag, K3 III) and HR6770 (V=4.64 
mag, G8 III). 
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Table 4 



Ground-based (M-D-M 2.4-m + Modspec) Observations 



Name 


Ubs. Date " 




Jixposure " 


Setup J 


NGC 


MM/YY 





(sec) 


Code 


— sm 


— I 1 /Qg; 


n — 


a x/ i Ann — 

X 10UU 





821 


11/95 


85 


5x1600 


6 


2778 


11/97 





23x1200 


5 


2778 


11/97 


70 


20x1200 


5 


3384 


03/97 





8x1200 


4 


3384 


01/97 


10 


6x1200 


6 


3384 


03/97 


70 


7x1200 


4 


3384 


01/97 


88 


6x1200 


6 


3608 


03/97 





7x1200 


4 


3608 


03/97 


75 


9x1200 


4 


4291 


03/97 





8x1200 


4 


4291 


03/97 


29 


2x1200 


4 


4291 


04/97 


70 


4x1200 


4 


4473 


03/97 





12x1200 


4 


4473 


03/97 


22 


6x1200 


4 


4473 


01/97 





4x1200 


6 


4473 


01/97 


90 


7x1200 


6 


4564 


03/98 





12x1200 


5 


4564 


03/98 


30 


8x1200 


5 


4564 


03/98 


70 


12x1200 


5 


4649 


01/97 





9x1200 


6 


4649 


03/97 


22 


8x1200 


4 


4649 


01/97 


28 


6x1200 


6 


4649 


01/97 


68 


9x1200 


6 


4697 


04/97 





9x1200 


4 


4697 


04/97 


48 


9x1200 


4 


4697 


04/97 


70 


18x1200 


4 


7457 


11/97 





20x1200 


5 


7457 


11/97 


70 


6x1200,6x1500 


4 



a The month and year of observation. 

b The position angle of the slit relative to the major axis of 
the galaxy. 

c Number of exposures x exposure length. 

d Code from Table 2 for the spectrograph setup used. 
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Table 5 



Nuker-law Profile fits 



NUC No. 


hilter " 






a 




'r 


Source t 


821 


H 


13.2 


— rrr 


1.0 




0.64 


Rav 


2778 

2778 


1 

V 


14.5 
15.1 


0.3 
0.1 


0.3 
0.4 


1.8 
1.6 


0.30 
0.16 


Nuk 
Nuk 


3384 
3384 


1 

V 


15.1 
16.3 


2.9 
2.9 


11.2 
12.0 


1.8 
1.8 


0.65 
0.64 


Nuk 
Nuk 


3608 


V 


15.5 


0.3 


1.1 


1.3 


0.00 


Fab 


4291 
4291 


H 

R 


12.3 
15.1 


0.5 
0.6 


2.1 
1.4 


1.5 
1.6 


0.02 
0.00 


Rav 
Rest 


4473 
4473 


1 

V 


16.2 
17.7 


5.3 
6.3 


0.7 
0.7 


2.7 
2.9 


0.01 
0.01 


Nuk 
Nuk 


4564 


R 


15.9 


1.3 


1.4 


1.3 


0.80 


Rest 


4649 


V 


17.2 


3.6 


2.0 


1.3 


0.15 


Fab 


4697 


V 


13.5 


2.6 


24.9 


1.0 


0.74 


Fab 


7457 


H 


13.5 


0.3 


2.3 


1.0 


0.35 


Rav 



a Filterused: H=F160W(nic), I=F814W, R=F702W, V=F555W. 
b Surface brightness at the break radius. 

c Break radius: the radius of the transition between inner and outer 
power-laws. 

d Inner power-law slope. 

e Source for fit parameters. Rav=Ravindranath et al. 2001; 
Fab=Faber et al. 1997; Rest=Rest et al. 2001; Nuk=Lauer 2002, pri- 
vate communication. 

Note. — These are fits to the surface brightness profile along the 
major axis. 



Table 6 



Bins for spectral extraction 



Name 


W") 


kange(pix) 


Width(pix) 


^ lib CCD 


cen 


0.000 


511 -511 


1 


rl 


0.051 


512-512 


1 


r2 


0.101 


513-513 


1 


r3 


0.177 


514-515 


2 


r4 


0.304 


516-518 


3 


r5 


0.507 


519-523 


5 


r6 


0.811 


524 - 530 


7 


11 


-0.051 


510-510 


1 


12 


-0.101 


509 - 509 


1 


13 


-0.177 


507 - 508 


2 


14 


-0.304 


504 - 506 


3 


15 


-0.507 


499 - 503 


5 


16 


-0.811 


492 - 498 


7 


Wilbur 


cen 


0.0 


300 - 300 


1 


rl 


0.37 


301 - 301 


1 


r2 


0.73 


302 - 302 


1 


r3 


1.28 


303 - 304 


2 


r4 


2.01 


305 - 306 


2 


r5 


3.10 


307- 310 


4 


r6 


4.93 


311 - 316 


6 


r7 


7.85 


317-325 


9 


r8 


12.23 


326 - 339 


14 


Charlotte 


cen 


0.0 


300 - 300 


1 


rl 


0.59 


301 - 301 


1 


r2 


1.48 


302 - 303 


2 


r3 


2.95 


304 - 306 


3 


r4 


4.72 


307 - 309 


3 


r5 


7.38 


310-315 


6 


r6 


11.80 


316-324 


9 


r7 


18.29 


326 - 336 


11 


r8 


28.32 


337 - 359 


23 


r9 


44.25 


360 - 390 


31 


rlO 


69.03 


391 -443 


52 



Name 
NGC 

82T 

2778 
3384 
3608 
4291 
4473 
4564 
4649 
4697 
7457 



Table 7 
Summary of STIS Kinematics 



&0,STIS 

km s -1 



248 (16) 
200 (29) 
128 ( 8) 
296 (27) 
306 (35) 
162 (25) 
220 (20) 
547(101) 
171 (10) 
105 ( 8) 



km s 
103 (23) 
61 (18) 
85 (10) 
28 (20) 
74(18) 
67(13) 
90 (34) 
145 (70) 
110(21) 
8(7) 



V grad 

km s _1 arcsec _1 



0.42 (0.10) 
0.43 (0.09) 
0.83 (0.08) 
0.09 (0.07) 
0.26 (0.04) 
0.46 (0.08) 
0.67 (0.06) 
0.27 (0.14) 
0.64 (0.13) 
0.38 (0.07) 



-5TTT 

129 
414 
113 
410 
278 
218 
403 
374 
11 



-0.04 (0.05 ) — -0.05 (0.03 ) 



-0.15(0.17) 
-0.04 (0.07) 
-0.13(0.08) 
-0.00 (0.07) 
-0.06 (0.04) 
0.00 (0.03) 
-0.07 (0.08) 
0.14(0.06) 
0.01 (0.05) 



-0.08 (0.08) 
-0.06 (0.03) 
-0.08 (0.04) 
-0.04 (0.05) 
-0.06 (0.03) 
-0.01 (0.03) 
-0.07 (0.04) 
0.10(0.03) 
-0.03 (0.03) 



a Dispersion measured in the central STIS bin with error from LOS VD fitting. 

b Maximum rotational velocity from symmetrized STIS data. Points with large scatter were ignored 
in the case of NGC 821 and NGC 2778. 

c Ratio of maximum rotation (over all data) and central STIS dispersion. 

d Velocity "gradient" is measured as the maximum of vl7/0.17 or v30/0.3, where vl7 is the sym- 
metrized rotation velocity at 0. 17", and v30 is the rotation at 0.30". The entire rotation curve is used 
to define the galaxy's systemic velocity, not the central bin. 

e The value of h4 in the central STIS bin (r=0."0). 

f The weighted mean of the 3 central h4 values (r=0"0, r=±0."05, not symmetrized). 
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Name 
NGC 

82T~ 

2778 
3384 
3608 
4291 
4473 
4564 
4649 
4697 
7457 



Table 8 

Summary of Modspec Kinematics 



&o,gb 
km s" 1 



a 

km s~ 



235 (3 
202(10. 
134(2) 
181 (5 
262(12 
183 (7) 
171 6 
373 ( 9) 
157(5' 
73(3 



192 (21) 
179 (52) 
119(7) 
150(12 
183(18 
197(10) 
112(17) 
318(37) 
145 (5) 
63 (3) 



C e C V max fiB V max /a 6 Wmaxfv ) 

km s" 1 km s" 1 



Q,CB 



<h'i> 1 



-70T 
175 
143 
182 
242 
190 
162 
375 
177 
67 



0.39 (0.05) 
0.47 (0.16) 
0.90 (0.09) 
0.16(0.04) 
0.43 (0.07) 
0.38 (0.03 
1.31 (0.20 
0.29 (0.09) 
0.70 (0.09' 
0.63 (0.12 



TX54- 
0.85 
0.99 
0.31 
0.71 
0.47 
1.18 
0.65 
0.84 
0.69 



""034 — -0.07 (0.02 



-0.03 (0.04) 
-0.04 (0.07) 
-0.04(0.10) 
-0.01 (0.03) 
-0.05 (0.02) 
-0.11 (0.02) 
-0.02 (0.04) 
-0.04 (0.02) 
-0.15(0.03) 
-0.02 (0.04) 



85 (14 
107 ( 8) 

24(5 
79 (10) 

75(4 
147 ( 4) 
93 (27) 
101 (13) 

40(7) 



0.24 
0.45 
0.21 
0.27 
0.40 
0.55 
0.16 
0.41 
0.46 



-0.06 (0.04 
-0.05 (0.01) 
-0.08 (0.02) 
-0.04 (0.03) 
-0.00 (0.02 
-0.01 (0.02 
0.00 (0.02) 
-0.04 (0.03' 
0.01 (0.03 



a Dispersion measured in central Modspec bin with error from LOS VD fitting. 

b Dispersion averaged along the Modspec slit, major axis, for r > 2."0, using symmetrized LOSVDs and an w l."0 slit 
width. The error is the standard deviation of averaged values. Only three numbers are averaged for NGC 2778 (the 6 
measurements outside of r =2."0 are symmetrized). The other galaxies include 4-8 numbers in the average. 

c Effective dispersion: the rms velocity relative to the systemic velocity of the galaxy, averaged over a slit width of ~ 
1" extending to R e . 

d Maximum rotational velocity from symmetrized Modspec data. Implausible maxima were discarded in the case of 
NGC 821, 4649, 4697. 

e Ratio of maximum rotation to a. The v max comes from the ground-based data. 

f The ratio of v max la to the v/W expected for an oblate, isotropic rotator calculated as ^/e/(l-e), where e is the isophote 
ellipticity. 

SEllipticity (l-b/a) from Bender (1988), Faber et al. (1997) (NGC 7457), Lauer et al. (1995) (NGC 3608), and Peletier 
et al. (1990) (NGC 2778). We measured e for NGC 4291 and 3384 using WFPC2 F555W images. 

h Fourth Gauss-Hermite coefficient from the r=0."0 bin of the major axis spectrum. 

'Third Gauss-Hermite coefficient averaged over r > 2"0 using the symmetrized, major axis LOSVD. The error is the 
standard deviation of the mean. 
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Table 9 

Ground-based (M-D-M 2.4-m + Modspec) Kinematics 



IN time 


PA / 


Sym7 b 


1 V l 


v( km s 1 ) J 


— 


a ( km s ! ) L 


■ — 


h3' 




h4 ^ 




— MS 9 1 f~\ftn\ 

INoZl^IVlgJ 




\es 


U.U 


7 

— / . 


71 

4. 


— 

Zjj. 


5 — 

j. 


m — 

u.uz 


— O 09 — 

U.UZ 


— 07 — 

—u.u / 


— O 09 — 

U.UZ 


iNoZl(IVlg) 


y 


var- 
ies 


U.D 


Zj. 


"7 


ZJ1. 


0. 


n fi9 

—U.UZ 


n ni 
U.Uj 


n m 
— U.U1 


n ni 
U.Uj 


INoZl^IVlg; 


X 


Vac 


1 9 

1 .Z 


^ 1 
j 1 . 


0. 


9 1 ^ 
Zl J. 


0. 


n 07 
—u.u / 


n n9 
U.UZ 


09 
U.UZ 


n (vx 
U.UJ 


W89 1 flVTaA 
IN OZ1 ^IVlgJ 




Vac 


1 .o 


S8 

JO. 


c 
j. 


90^ 
ZUj. 


7 


n 07 
u.u / 


CiX 

U.UJ 


OS 
U.UJ 


03 
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9. 
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9. 
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0.02 


0.00 


0.02 


N4649 


22 


Yes 


0.6 


12. 


8. 


382. 


9. 


-0.01 


0.02 


-0.02 


0.02 


N4649 


22 


Yes 


1.5 


9. 


8. 


371. 


10. 


-0.02 


0.02 


0.01 


0.02 


N4649 


22 


Yes 


3.0 


20. 


8. 


360. 


10. 


-0.05 


0.02 


0.01 


0.02 


N4649 


22 


Yes 


4.7 


15. 


7. 


355. 


11. 


-0.05 


0.02 


0.02 


0.03 


N4649 


22 


Yes 


7.4 


20. 


7. 


323. 


10. 


-0.04 


0.02 


-0.03 


0.02 


N4649 


22 


Yes 


11.8 


43. 


11. 


291. 


15. 


-0.03 


0.03 


0.04 


0.03 


N4649 


22 


Yes 


18.3 


69. 


17. 


280. 


22. 


-0.03 


0.05 


0.04 


0.04 


N4649 


22 


Yes 


28.3 


103. 


38. 


285. 


69. 


-0.02 


0.09 


0.04 


0.10 
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28 
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68 


Yes 
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4. 


8. 


353. 


7. 


0.00 


0.02 


-0.11 


0.01 


N4649(Mg) 


68 


Yes 


3.0 


10. 


8. 


341. 


8. 


-0.01 


0.02 


-0.09 


0.01 


N4649(Mg) 


68 


Yes 


4.7 


-1. 


8. 


330. 


6. 


-0.04 


0.02 


-0.03 


0.01 


N4649(Mg) 


68 


Yes 


7.4 


8. 


5. 


327. 


7. 


-0.02 


0.01 


-0.01 


0.01 


N4649(Mg) 


68 


Yes 


11.8 


14. 


6. 


316. 


6. 


0.01 


0.01 


0.01 


0.01 


N4649(Mg) 


68 


Yes 


18.3 


43. 


8. 


294. 


6. 


0.02 


0.03 


-0.02 


0.01 


N4649(Mg) 


68 


Yes 


28.3 


36. 


10. 


291. 


9. 


-0.01 


0.03 


-0.02 


0.02 


N4649(Mg) 


68 


Yes 


44.2 


52. 


10. 


239. 


10. 


0.05 


0.03 


-0.03 


0.02 


N4697 





Yes 


0.0 


-22. 


5. 


157. 


5. 


0.02 


0.03 


-0.04 


0.03 


N4697 





Yes 


0.6 


28. 


3. 


163. 


4. 
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0.02 


-0.01 


0.02 



N4697 Yes 1.2 61. 3. 149. 4. -0.08 0.02 0.02 0.02 

xt i /ir\i f\ ^/-_„ ^ 1 in o i i n a n 1 1 r\ f\n n ni n nn 



41 



TABLE 9 — Continued 



Name 




Sym'.' b 


r{") L 


v(kms ') J 




<r(kms l f 




h3' 




h4B 


e*4 


N7457 





Yes 


2.0 


=T9— 




64. 


5. 


0.01 


0.04 


-0.02 


0.03 


N7457 





Yes 


3.1 


-26. 


4. 


58. 


3. 


0.00 


0.03 


-0.06 


0.02 


N7457 





Yes 


4.9 


-33. 


6. 


63. 


7. 


-0.08 


0.04 


-0.02 


0.03 


N7457 





Yes 


7.7 


-40. 


7. 


65. 


7. 


-0.01 


0.05 


-0.04 


0.02 


N7457 





No 


0.00 


1. 


3. 


73. 


3. 


-0.05 


0.03 


0.01 


0.03 


N7457 





No 


-0.37 


3. 


3. 


67. 


2. 


0.00 


0.03 


-0.06 


0.01 


N7457 





No 


-0.73 


4. 


4. 


67. 


4. 


0.00 


0.04 


-0.03 


0.02 


N7457 





No 


-1.28 


16. 


3. 


82. 


4. 


-0.01 


0.04 


0.04 


0.05 


N7457 





No 


-2.01 


12. 


5. 


70. 


6. 


0.00 


0.06 


-0.05 


0.02 


N7457 





No 


-3.10 


24. 


6. 


47. 


7. 


-0.01 


0.04 


-0.07 


0.01 


N7457 





No 


-4.93 


32. 


8. 


85. 


10. 


0.17 


0.08 


0.06 


0.07 


N7457 





No 


-7.85 


52. 


10. 


59. 


9. 


-0.02 


0.05 


-0.05 


0.03 


N7457 





No 


0.37 


-3. 


3. 


74. 


3. 


-0.03 


0.04 


-0.03 


0.02 


N7457 





No 


0.73 


-7. 


4. 


68. 


4. 


-0.06 


0.04 


0.02 


0.03 


N7457 





No 


1.28 


-15. 


4. 


61. 


5. 


0.02 


0.05 


-0.01 


0.04 


N7457 





No 


2.01 


-21. 


5. 


55. 


8. 


0.00 


0.05 


-0.06 


0.02 


N7457 





No 


3.10 


-24. 


6. 


74. 


8. 


0.09 


0.08 


0.01 


0.06 


N7457 





No 


4.93 


-40. 


8. 


60. 


5. 


-0.02 


0.04 


-0.05 


0.02 


N7457 





No 


7.66 


-28. 


10. 


77. 


10. 


-0.06 


0.06 


-0.05 


0.04 


N7457 


70 


Yes 


0.0 


-2. 


4. 


65. 


3. 


-0.09 


0.04 


0.01 


0.03 


N7457 


70 


Yes 


0.4 


2. 


2. 


51. 


4. 


-0.08 


0.04 


-0.04 


0.01 


N7457 


70 


Yes 


0.7 


0. 


3. 


57. 


5. 


-0.02 


0.05 


-0.05 


0.02 


N7457 


70 


Yes 


1.3 


1. 


4. 


68. 


5. 


-0.04 


0.05 


-0.01 


0.02 


N7457 


70 


Yes 


2.0 


-4. 


6. 


60. 


5. 


-0.04 


0.04 


-0.03 


0.02 


N7457 


70 


Yes 


3.1 


-3. 


8. 


57. 


6. 


-0.01 


0.05 


-0.02 


0.02 



"Position angle of slit relative to the major axis. 

b If "Yes", the spectra from the two sides of the galaxy center have been symmetrized. 
c Radius, in arcseconds. 

d Line-of-sight velocity, in km s _1 , from fit to Equation 1 (see text). 
e Velocity dispersion, km s -1 , from fit to Equation 1 (see text). 
f Third Gauss-Hermite coefficient. 
E Fourth Gauss-Hermite coefficient. 



Table 10 

stis-measured kinematics from ca ii 
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n ri+T 
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hi 1 
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h4 ^ 


e h4 


N821 
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N 


— 0TO0 


-1. 
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238— 
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0.02 


0.05 


— 0.04 


0.05 


N821 
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0.05 
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0.06 


0.05 
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0.05 


N821 
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0.10 


-101. 


20. 


244. 
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0.05 


0.05 


-0.02 


0.07 


N821 


-150. 
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0.17 
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186. 
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0.04 


0.03 


0.00 


0.05 


N821 
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N 


0.30 
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24. 


177. 
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0.10 


0.07 


0.02 


0.08 


N821 
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0.50 


-47. 


32. 


206. 
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0.05 


0.07 
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0.08 


N821 


-150. 
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0.82 
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233. 
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0.09 


0.15 


0.04 


0.16 


N821 
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N 


-0.05 


5. 


20. 


229. 


16. 


-0.08 


0.05 


-0.05 


0.05 


N821 


-150. 


N 


-0.10 


62. 


20. 


218. 


18. 


-0.11 


0.06 


-0.08 


0.05 


N821 


-150. 


N 


-0.17 


68. 


25. 


235. 


20. 


-0.14 


0.06 


-0.10 


0.05 


N821 


-150. 


N 


-0.30 


194. 


30. 


222. 


23. 


-0.19 


0.09 


0.10 


0.10 


N821 


-150. 


N 


-0.50 


116. 


47. 


191. 


33. 


0.04 


0.06 


-0.06 


0.06 


N821 


-150. 


N 


-0.82 


231. 


45. 


238. 


24. 


-0.23 


0.10 


0.00 


0.12 


N821 


-150 


Y 


0.00 


-1. 


19. 


248. 


16. 


0.02 


0.05 


-0.04 


0.05 


N821 


-150 


Y 


0.05 


33. 


12. 


244. 


12. 


-0.04 


0.04 


-0.09 


0.04 


N821 


-150 


Y 


0.10 


83. 


14. 


217. 


15. 


-0.11 


0.04 


-0.09 


0.04 


N821 


-150 


Y 


0.17 


87. 


16. 


192. 


13. 


-0.14 


0.03 


0.01 


0.04 


N821 


-150 


Y 


0.30 


103. 


23. 


209. 


25. 


-0.11 


0.06 


0.06 


0.06 


N821 


-150 


Y 


0.50 


65. 


28. 


197. 


29. 


-0.02 


0.06 


-0.06 


0.04 


N821 


-150 


Y 


0.82 


189. 


33. 


228. 


24. 


-0.16 


0.09 


-0.04 


0.08 


N2778 


50 


N 


0.00 


-66. 


39. 


200. 


29. 


0.09 


0.13 


-0.15 


0.17 


N2778 


50 


N 


0.05 


32. 


34. 


149. 


30. 


-0.03 


0.07 


-0.07 


0.10 


N2778 


50 


N 


0.10 


67. 


31. 


158. 


34. 


-0.05 


0.12 


-0.11 


0.11 


N2778 


50 


N 


0.17 


-35. 


45. 


194. 


38. 


-0.10 


0.14 


-0.05 


0.16 


N2778 


50 


N 


0.43 


54. 


26. 


99. 


33. 


-0.02 


0.05 


-0.08 


0.06 


N2778 


50 


N 


-0.05 


-152. 


71. 


206. 


64. 


-0.03 


0.19 


-0.03 


0.16 


N2778 


50 


N 


-0.10 


-63. 


37. 


153. 


26. 


-0.08 


0.09 


-0.01 


0.08 


N2778 


50 


N 


-0.17 


-48. 


38. 


128. 


25. 


0.05 


0.10 


-0.03 


0.07 


N2778 


50 


N 


-0.43 


-60. 


36. 


126. 


27. 


0.13 


0.12 


0.00 


0.13 


N2778 


50 


Y 


0.00 


-66. 


39. 


200. 


29. 


0.09 


0.13 


-0.15 


0.17 


N2778 


50 


Y 


0.05 


-45. 


25. 


183. 


38. 


0.14 


0.10 


0.00 


0.12 


N2778 


50 


Y 


0.10 


-61. 


18. 


142. 


17. 


0.00 


0.07 


-0.11 


0.05 


N2778 


50 


Y 


0.17 


-15. 


23. 


155. 


24. 


0.08 


0.09 


-0.10 


0.09 


N2778 


50 


Y 


0.43 


-55. 


19. 


113. 


19. 


0.07 


0.08 


-0.07 


0.06 


N3384 


58 


N 


0.00 


-15. 


9. 


128. 


8. 


0.01 


0.06 


-0.04 


0.07 


N3384 


58 


N 


-0.05 


16. 


8. 


133. 


8. 


-0.06 


0.06 


-0.10 


0.06 


N3384 


58 


N 


-0.10 


48. 


7. 


109. 


9. 


-0.14 


0.08 


0.09 


0.06 


N3384 


58 


N 


-0.17 


79. 


6. 


84. 


18. 


-0.11 


0.07 


0.05 


0.11 


N3384 


58 


N 


-0.30 


75. 


10. 


98. 


16. 


-0.12 


0.09 


-0.01 


0.11 


N3384 


58 


N 


-0.50 


87. 


10. 


99. 


14. 


-0.13 


0.09 


0.01 


0.08 


N3384 


58 


N 


-0.80 


95. 


13. 


95. 


21. 


-0.01 


0.06 


-0.04 


0.06 


N3384 


58 


N 


0.05 


-53. 


7. 


122. 


6. 


0.02 


0.05 


-0.04 


0.05 


N3384 


58 


N 


0.10 


-66. 


9. 


104. 


10. 


-0.04 


0.06 


0.03 


0.05 


N3384 


58 


N 


0.17 


-65. 


9. 


86. 


9. 


0.07 


0.07 


-0.03 


0.06 


N3384 


58 


N 


0.30 


-62. 


19. 


121. 


22. 


-0.06 


0.11 


0.06 


0.12 


N3384 


58 


N 


0.50 


-54. 


16. 


107. 


15. 


0.08 


0.08 


0.00 


0.07 


N3384 


58 


N 


0.80 


-72. 


22. 


98. 


17. 


0.04 


0.09 


-0.06 


0.06 


N3384 


58 


Y 


0.00 


-15. 


9. 


128. 


8. 


0.01 


0.06 


-0.04 


0.07 


N3384 


58 


Y 


0.05 


35. 


5. 


127. 


4. 


-0.05 


0.04 


-0.08 


0.04 


N3384 


58 


Y 


0.10 


57. 


7. 


102. 


7. 


-0.04 


0.05 


0.07 


0.06 


N3384 


58 


Y 


0.17 


70. 


5. 


85. 


4. 


-0.10 


0.03 


-0.02 


0.03 


N3384 


58 


Y 


0.30 


69. 


9. 


108. 


10. 


-0.06 


0.08 


0.00 


0.09 


N3384 


58 


Y 


0.50 


74. 


10. 


101. 


9. 


-0.15 


0.07 


0.05 


0.07 


N3384 


58 


Y 


0.80 


85. 


10. 


92. 


11. 


-0.04 


0.06 


-0.06 


0.06 


N36U8 


94 


N 


0.00 


11. 


30. 


296. 


27. 


0.05 


0.08 


-0.13 


0.08 


N3608 


94 


N 


-0.05 


28. 


23. 


270. 


25. 


-0.02 


0.07 


-0.04 


0.08 


N3608 


94 


N 


-0.10 


-1. 


23. 


266. 


21. 


0.03 


0.07 


-0.10 


0.07 


N3608 


94 


N 


-0.17 


-30. 


26. 


289. 


35. 


-0.02 


0.10 


0.04 


0.08 


N3608 


94 


N 


-0.30 


10. 


27. 


227. 


28. 


0.10 


0.08 


-0.01 


0.08 


N3608 


94 


N 


-0.50 


5. 


20. 


193. 


18. 


0.05 


0.06 


0.01 


0.04 


N3608 


94 


N 


-0.80 


51. 


43. 


168. 


38. 


0.01 


0.11 


-0.04 


0.14 


N3608 


94 


N 


0.05 


17. 


27. 


256. 


24. 


0.05 


0.07 


-0.08 


0.07 


N3608 


94 


N 


0.10 


-14. 


25. 


206. 


21. 


-0.02 


0.07 


-0.07 


0.06 


N3608 


94 


N 


0.17 


8. 


25. 


265. 


38. 


0.03 


0.11 


0.08 


0.10 


N3608 


94 


N 


0.30 


-11. 


23. 


187. 


21. 


0.05 


0.07 


-0.06 


0.07 


N3608 


94 


N 


0.50 


-7. 


27. 


238. 


33. 


0.19 


0.11 


0.14 


0.10 


N3608 


94 


N 


0.80 


2. 


36. 


161. 


38. 


0.05 


0.10 


0.00 


0.13 


N3608 


94 


Y 


0.00 


11. 


30. 


296. 


27. 


0.05 


0.08 


-0.13 


0.08 


N3608 


94 


Y 


0.05 


2. 


14. 


256. 


15. 


-0.02 


0.05 


-0.08 


0.03 


N3608 


94 


Y 


0.10 


11. 


15. 


233. 


15. 


0.02 


0.05 


-0.09 


0.03 


N3608 


94 


Y 


0.17 


-19. 


19. 


271. 


31. 


-0.03 


0.07 


0.07 


0.05 


N3608 


94 


Y 


0.30 


5. 


14. 


194. 


15. 


-0.02 


0.05 


-0.06 


0.03 


N3608 


94 


Y 


0.50 


14. 


13. 


192. 


17. 


-0.03 


0.05 


0.07 


0.04 


N3608 


94 


Y 


0.80 


28. 


20. 


148. 


27. 


0.01 


0.06 


-0.03 


0.05 


N4291 


110 


N 


0.00 


51. 


29. 


306. 


35. 


0.05 


0.08 


0.00 


0.07 


N4291 


110 


N 


0.05 


-76. 


40. 


310. 


38. 


0.09 


0.09 


-0.01 


0.08 


N4291 


110 


N 


0.10 


-118. 


38. 


304. 


38. 


0.15 


0.10 


0.05 


0.11 


N4291 


110 


N 


0.17 


-104. 


24. 


264. 


29. 


-0.02 


0.07 


0.03 


0.06 


N4291 


110 


N 


0.30 


-114. 


32. 


268. 


33. 


-0.07 


0.08 


-0.03 


0.06 


N4291 


110 


N 


0.50 


-48. 


27. 


260. 


18. 


0.04 


0.06 


-0.08 


0.04 


N4291 


110 


N 


0.80 


-74. 


32. 


279. 


33. 


0.12 


0.10 


0.07 


0.07 


N4291 


110 


N 


-0.05 


47. 


33. 


354. 


45. 


0.05 


0.10 


0.02 


0.08 


N4291 


110 


N 


-0.10 


-47. 


25. 


261. 


34. 


0.06 


0.09 


0.01 


0.08 


N4291 


110 


N 


-0.17 


34. 


21. 


275. 


21. 


-0.01 


0.06 


-0.06 


0.06 


N4291 


110 


N 


-0.30 


15. 


27. 


273. 


22. 


-0.07 


0.07 


-0.05 


0.06 


N4291 


110 


N 


-0.50 


97. 


26. 


253. 


19. 


0.00 


0.04 


-0.08 


0.04 


N4291 


110 


N 


-0.80 


10. 


34. 


278. 


34. 


-0.02 


0.07 


-0.06 


0.05 


N4291 


110 


Y 


0.00 


51. 


29. 


306. 


35. 


0.05 


0.08 


0.00 


0.07 


N4291 


110 


Y 


0.05 


49. 


22. 


282. 


20. 


-0.11 


0.06 


-0.07 


0.06 


N4291 


110 


Y 


0.10 


19. 


21. 


262. 


20. 


-0.05 


0.05 


-0.07 


0.04 


N4291 


110 


Y 


0.17 


70. 


18. 


271. 


20. 


0.00 


0.06 


0.01 


0.05 


N4291 


110 


Y 


0.30 


63. 


18. 


252. 


16. 


-0.02 


0.04 


-0.05 


0.03 


N4291 


110 


Y 


0.50 


74. 


18. 


241. 


16. 


-0.04 


0.05 


-0.11 


0.03 


N4291 


110 


Y 


0.80 


61. 


23. 


246. 


24. 


-0.06 


0.06 


-0.03 


0.06 


N4473 


-100 


N 


0.00 


-A. 


25. 


162. 


25. 


-0.02 


0.07 


-0.08 


0.08 


N4473 


-100 


N 


-0.05 


-A. 


17. 


153. 


14. 


0.05 


0.05 


-0.05 


0.05 


N4473 


-100 


N 


-0.10 


27. 


14. 


140. 


12. 


-0.05 


0.04 


0.01 


0.04 


N4473 


-100 


N 


-0.17 


46. 


17. 


150. 


16. 


-0.04 


0.05 


-0.03 


0.04 


N4473 


-100 


N 


-0.30 


69. 


19. 


155. 


18. 


-0.06 


0.05 


-0.02 


0.05 


N4473 


-100 


N 


-0.50 


67. 


20. 


157. 


13. 


-0.08 


0.04 


-0.04 


0.04 


N4473 


-100 


N 


-0.80 


57. 


17. 


199. 


15. 


-0.05 


0.05 


-0.03 


0.04 


N4473 


-100 


N 


0.05 


-43. 


20. 


163. 


14. 


0.03 


0.06 


-0.06 


0.04 


N4473 


-100 


N 


0.10 


-41. 


20. 


171. 


15. 


0.08 


0.06 


-0.06 


0.05 
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TABLE 10 — Continued 



Name 


PA" 


Sym'. ,lj 


r("f 


v(kms ') u 


e v 


a ( km s 1 r 


e CT 






h4 K 


e*4 


N4649 


-95 — 


Y 


1730- 


=T3— 




439— 




-0.05 


0.06 


0.01 


0.05 


N4697 


68 


N 


0.00 


-7. 


10. 


171. 


10. 


-0.05 


0.04 


0.10 


0.04 


N4697 


68 


N 


0.05 


-28. 


10. 


177. 


11. 


0.04 


0.04 


0.09 


0.05 


N4697 


68 


N 


0.10 


-21. 


12. 


215. 


10. 


0.06 


0.04 


-0.02 


0.05 


N4697 


68 


N 


0.17 


-54. 


16. 


206. 


22. 


0.05 


0.07 


0.05 


0.08 


N4697 


68 


N 


0.30 


-48. 


16. 


153. 


36. 


-0.01 


0.12 


0.02 


0.15 


N4697 


68 


N 


0.50 


-63. 


32. 


149. 


41. 


-0.15 


0.14 


0.11 


0.17 


N4697 


68 


N 


0.80 


-104. 


24. 


115. 


28. 


-0.05 


0.10 


-0.06 


0.11 


N4697 


68 


N 


-0.05 


9. 


8. 


168. 


11. 


-0.02 


0.04 


0.10 


0.04 


N4697 


68 


N 


-0.10 


28. 


11. 


182. 


15. 


0.07 


0.04 


0.07 


0.05 


N4697 


68 


N 


-0.17 


77. 


16. 


171. 


20. 


0.10 


0.08 


0.09 


0.07 


N4697 


68 


N 


-0.30 


54. 


16. 


166. 


21. 


0.03 


0.11 


-0.03 


0.07 


N4697 


68 


N 


-0.50 


82. 


27. 


115. 


37. 


0.22 


0.16 


0.16 


0.16 


N4697 


68 


N 


-0.80 


155. 


41. 


201. 


38. 


0.46 


0.14 


0.28 


0.18 


N4697 


68 


Y 


0.00 


0. 


0. 


155. 


13. 


0.00 


0.00 


0.14 


0.06 


N4697 


68 


Y 


0.05 


20. 


7. 


187. 


20. 


-0.02 


0.02 


0.16 


0.07 


N4697 


68 


Y 


0.10 


27. 


7. 


217. 


9. 


0.00 


0.03 


0.04 


0.06 


N4697 


68 


Y 


0.17 


64. 


10. 


195. 


18. 


0.00 


0.04 


0.11 


0.08 


N4697 


68 


Y 


0.30 


51. 


14. 


178. 


20. 


0.02 


0.08 


0.04 


0.09 


N4697 


68 


Y 


0.50 


67. 


25. 


126. 


39. 


0.19 


0.16 


0.27 


0.12 


N4697 


68 


Y 


0.80 


110. 


21. 


138. 


24. 


0.16 


0.16 


-0.06 


0.17 


N7457 


130 


N 


0.00 


-3. 


9. 


105. 


8. 


-0.07 


0.06 


0.01 


0.05 


N7457 


130 


N 


0.05 


-8. 


9. 


100. 


8. 


-0.05 


0.07 


-0.06 


0.04 


N7457 


130 


N 


0.10 


8. 


17. 


82. 


14. 


-0.11 


0.07 


-0.04 


0.10 


N7457 


130 


N 


0.17 


-19. 


41. 


92. 


23. 


-0.09 


0.15 


-0.11 


0.22 


N7457 


130 


N 


-0.05 


9. 


11. 


79. 


7. 


-0.08 


0.05 


-0.03 


0.04 


N7457 


130 


N 


-0.10 


-11. 


12. 


76. 


8. 


-0.15 


0.09 


0.05 


0.06 


N7457 


130 


N 


-0.17 


-10. 


19. 


77. 


13. 


-0.11 


0.07 


-0.01 


0.06 


N7457 


130 


Y 


0.00 


-3. 


9. 


105. 


8. 


-0.07 


0.06 


0.01 


0.05 


N7457 


130 


Y 


0.05 


8. 


7. 


88. 


6. 


-0.05 


0.05 


-0.03 


0.03 


N7457 


130 


Y 


0.10 


-7. 


7. 


80. 


10. 


-0.01 


0.06 


-0.03 


0.05 


N7457 


130 


Y 


0.17 


-2. 


15. 


79. 


18. 


0.06 


0.08 


-0.03 


0.07 



"Position angle of STIS slit in degrees. 

b If "Y", the spectra on the two sides of the galaxy center have been symmetrized. 
c Radius, in arcseconds. 

d Line-of-sight velocity, in km s -1 , from fit to Equation 1 . 
e Velocity dispersion, km s~' , from fit to Equation 1. 
'Third Gauss-Hermite coefficient. 
g Fourth Gauss-Hermite coefficient. 



